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NOTICE 

This report was prepared as an account of Government sponsored 
work. Neither the United States, nor the National Aeronautics 
and Space Administration (NASA), nor any person acting on be- 
half of NASA: 

A. ) Makes any warranty or representation, expressed or 

implied, with respect to the accuracy, completeness, 
or usefulness of the information contained in this report, 
or that the use of any information, apparatus, method, 
or process disclosed in this report may not infringe 
privately owned rights; or 

B. ) Assumes any liabilities with respect to the use of, or 

for damages resulting from the use of any information, 
apparatus, method or process disclosed in this report. 

As used above, "person acting on behalf of NASA" includes any 
employee or contractor of NASA, or employee of such con- 
tractor, to the extent that such employee or contractor of NASA, 
or employee of such contractor prepares, disseminates, or provides 
access to, any information pursuant to his employment or contract 
with NASA, or his employment with such contractor. 
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1 .0 SUMMARY 

flu TH ~*9&$ 

This report covers the work performed between July 1 and November 30, 1962 on Contract 
NAS 3-2159, Wetting and Non-Wetting Mercury Condensing Research . The investigation 
consists of both an analytical and an experimental phase. The experimentation consists of 
local pressure drop measurements for condensing mercury in horizontal straight and tapered 
steel tubes. Both wetting and non-wetting conditions are explored. The analysis consists 
of the formulation of a dropwise condensing fluid mechanic model for correlating the data 
obtained from the non-wetting forced convection condensation tests. Also included is a 
hydrodynamic stability analysis of liquid films for application to wetting condensing and its 
multitude of possible flow patterns. This analysis is required to match the proper flow regime 
with its pressure drop correlation. 

The test data are compared with the Lockhart-Martinelli correlation to indicate agreement 
and note shortcomings. The state of the art of two-phase fluid mechanics is reviewed to 
give the reader a proper perspective of the work being done and the work that has been done. 

Non-wetting pressure drop data are tabulated herein with only preliminary attempts of 
correlation indicated. The final report will include the wetting data currently being 
obtained as well as the comprehensive results of correlation currently being compiled. 
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2.0 INTRODUCTION 


The purpose of the wetting and non-wetting condensing research program, sponsored by 
NASA under contract NAS3-2159, was to refine procedures for designing low heat flux* 
mercury condensers for space power systems. This study placed particular emphasis on the 
differences between wetting and non-wetting condenser performance and the different de- 
sign considerations required. 

Mercury is considered to be "non-wetting" (i.e. , the contact angle is greater than 90 degrees) 
with most fabrication materials presently planned for use in Rankine cycle space power plants. 
These systems will operate up to temperatures of approximately 1200°F or less (Ref 1 and 2), 
and presenrly are constructed from either the stainless steels or Haynes-25 type (high cobalt 
content) materials. It has been observed, however, that mercury will wet these materials 
with sufficient time and temperature. 

The change from non-wetting to wetting significantly affects the condensing flow regime, 
and therefore wetting and non-wetting must be understood by designers of mercury condensers. 
During non-wetting a dropwise flow regime has been obtained in glass tubes in tests run at 
TRW. Some over-all condensing and adiabatic data (e.g. , total pressure drop) have also been 
obtained at TRW. References 3, 4, and 5 present an analysis of dropwise condensation of 
mercury. Further development of this drop analysis was required to account for the change 
in drop size with length. In addition, local pressure measurements along the tube in more 
extensive temperature and flow ranges were required. 

When mercury wets the wall a thin film of liquid mercury forms on the tube wall. This film 
is more unstable than drops; consequently, more flow regimes exist and the fluid dynamics 
are more complex. When mercury forms a continous liquid film which can readily transmit 
wave phenomena, sensitivity to small body force disturbances may increase. The waves may 
cause spray (drop formation) as well as serious slugging with attendant pressure and inventory 
fluctuations. A significant amount of experimentation and some analysis has already been 
done on adiabatic wetting two-phase flow and usefulness to the case at hand was investigated. 
However, data with heat transfer, more analysis, and local pressure measurements were still 
required for the case of wetting mercury condensation. 

The program described in this report was instigated to obtain local pressure drop data for the 
wetting and non-wetting condensing of mercury. The state of the art of wetting and non- 
wetting condensing theory and analysis is established and refinements are made. 


* In direct radiator-condenser systems ^condensing heat flux is limited by the heat rejection 
area to approximately 15,000 Btu/hr-ft . 
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3.0 PROGRAM OBJECTIVES 


The basic program objectives were to provide both experimental data and mathematical 
tools for the design of a low-heat flux zero gravity mercury condenser In which either 
wetting or nonwetting condensing may occur. The data should be useful to check on 
theory as well as to provide raw data for the designer of condenser components. Both the 
data and analysts should be applicable to condenser materials, configurations, flow ranges, 
and saturation levels of Interest for present-day mercury Ranklne cycle space power systems. 

Stainless steel was originally selected as the condensing tube material, but Haynes-25 
type material was subsequently substituted. The geometries selected for test are listed 
below: 


Conditions Test 


Configuration 

Non-Wettlng 

Wetting 

Series 

0.319 inch ID constant diameter 
by 8 feet long 

Yes 

Yes 

A,F 

0.397 inch ID constant diameter 
by 8 feet long 

Yes 

No 

E 

0.4 to 0.2 inch ID linear 

Yes 

Yes 

D 


diameter taper by 7 feet long 


These geometries were chosen because they cover the range of tube sizes presently under 
consideration in advanced mercury direct radiator-condensers. 

The ranges of flow rate, pressure level, and condensing lengths over which data are of in- 
terest for direct application to advanced mercury condenser design are summarized below: 



Flow Rate 

Inlet Pressure 

Condensing 

Tube Configuration 

lb/min 

Psia 

Length, ft 

Constant Diameter 

1 1/4 to 3 

11 to 30 

4 to 8 

Tapered 

1 1/4 to 3 

11 to 30 

7.5 


A greater variation of flow rate, pressure level, condensing length and configuration Is 
also of Interest, as extreme variations of Reynold's and Froude numbers would be useful for 
extrapolation of the data. From system considerations, however, pressure levels higher 
than 30 psia have been found impractical as cycle efficiency is unnecessarily decreased. 
Also, for the tube sizes and saturation levels chosen, flow rate must be limited because 
unreasonably high vapor inlet velocities may result. Also,flows less than about 1 1/4 
pounds per minute for the tube sizes investigated are too sensitive to gravity to provide any 
useful zero-gravity data. To provide greater insight into the phenomena occurring. 
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pressure readings were taken at intervals of 18 inches along the condensing section. 
Pressure readings at smaller Intervals might affect the Internal phenomena while larger 
Intervals would be less "local". 

To establish and refine the present state of the art of condensing theory and analysis (1) 
existing adiabatic two-phase flow theory was applied to the case of wetting condensation; 
(2) film buildup and film stability during wetting condensation were considered; and (3) 
the existing theory for non-wetting pressure drop was refined to include the effect of drop 
size variation along the condensing tube. 
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4.0 TWO-PHASE FLOW THEORY AND ANALYSIS 


Design of a forced convection condenser for operation in a variable gravity environment 
requires consideration of the following factors.’ 

f. flow regime, 

2. pressure drop, 

3. density, or liquid and vapor contents 

4. amplitude and frequency of pressure fluctuation. 

Knowledge of the flow regime is required so that the proper pressure drop and density 
equations may be applied. Pressure drop must be calculated to insure that system require- 
ments are met. Inventory allocation must also be known and therefore mean density or 
liquid and vapor contents must be known. Finally, pressure fluctuations must be investi- 
gated to insure that excessive oscillations are not present. Tolerance limits for pressure 
fluctuations must first be defined by system (principally pump and turbine) requirements. 

The present state of the art of condenser design allows only estimation of the relationships 
between the above parameters on the basis of correlations of adiabatic test results, and the 
estimation is primarily true for wetting fluids only. The sections below review modified 
adiabatic approaches and present unique analytical treatments for the prediction of wetting 
film stability and non-wetting mercury condensing pressure drop. 

The design procedure for a zero-gravity condenser should be similiar to that for the horizontal 
1 g case, if the annular flow regime is predominant during wetting condensing (which it Is, 
as indicated in Appendix B) and if the drop agglomeration is small during the non-wetting. 
Horizontal 1 g condensing data in constant diameter tubes with high vapor inlet velocities 
and in tapered tubes with high velocities throughout are expected to be very similiar to 
zero-gravity data. The high axial velocity component creates an initial inertia force that 
predominates over the gravity force, thereby minimizing the gravity effect. Based on the 
comparison of some data with theory in Section 7.0, inlet vapor velocities of roughly 
150 feet per second or higher appear adequate. 

In a constant diameter tube, the major frictional and droplet or wave drag effects on pressure 
drop occur in the high vapor velocity or high quality region and therefore should not be 
significantly affected by gravity. In the low quality (velocity) region, however, liquid 
distribution is significantly affected by gravity. As a result, the pressure pulsations and 
local pressure drop may be significantly affected. The over-all pressure drop for the tube 
is not significantly affected, however, since the low vapor velocity should result in low 
pressure gradients. 

In tapered tubes the velocity is almost entirely constant. If the velocity is sufficiently 
high, identical flow regimes will result in the horizontal 1 g and zero gravity cases. The 
exact effects of gravity on the condensing process are only a matter of speculation at this 
time and analysis must await the completion of zero-gravity experimentation. 
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It has been observed that mercury systems can operate with either wetting or non-wetting 
condensation. The surface chemistry aspects relating the wetting and non-wetting states 
are complex and the transition is not well understood. Consequently, the wetting and non- 
wetting analytical considerations that are required are presented independently below. 

4. 1 CONDENSATION OF WETTING FLUIDS 

Two approaches to the prediction of the flow regime, pressure drop, density, and pressure 
fluctuations are possible for the wetting condensation of fluids such as mercury: 

1. The annular flow pattern can be assumed. Local filmthickness along the tube 
can be computed and energy dissipation and stability can be analyzed. A key 
missing link in this approach is the analytical tools required to perform the 
stability calculations. Section 4.1.1 reviews the situation. 

2. Existing two-phase adiabatic and diabatic data and correlations may be employed. 
Section 4.1.2, “Horizontal Two-Phase Flow," reviews this approach. 

The film stability approach described in Section 4.1.1 is more theoretically correct, but 
also is less understood and thus cannot be used for actual design of a wetting condenser at 
this time. Practically speaking, the "modified adiabatic" approach to the fluid dynamic design 
of a wetting mercury condenser suggested in Section 4.1.2 is the only means presently avail- 
able. The effects of heat transfer, zero gravity, and equilibrium length requirements in 
this approach are still unknown. 

The following sections contain much discussion of adiabatic two-phase phenomena, be- 
cause observing diabatic phenomena as precisely as required is difficult. Wherever possible, 
however, the considerations required to relate adiabatic and diabatic results are presented. 

4.1.1 Film Stability Analysis 

The information required to approach the design of a wetting condenser by analyzing film 
buildup and stability is far from complete. This is basically due to lack of (1) consistency 
and experimental verification of the theories for the stability of thin films and (2) means 
for predicting the effect of surface irregularities (waves, etc.) on pressure drop. In the 
section below, reason (1) is discussed extensively as the calculation of film stability is 
important in determining the flow regime, pressure drop, density, and amplitude and fre- 
quency of fluctuation. 

Pressure drop calculations for the case of an annular film could be made on the basis of the 
film wave friction factors for the wave heights and lengths that may exist. Very little in- 
formation of this sort is known to be available. In case liquid droplets are entrained in the 
vapor as a result of film breakup, a consideration of the drag effect of these drops in transit 
is required. The analysis of Section 4.2 would then be applicable. 

Computation of the local film thickness at each point along the tube requires determination 
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of three factors; 


1. Interfacial shear stress caused by condensing momentum transfer and frictional 
shear, 

2. body force vector, 

3. liquid flow rate. 

The local vapor and liquid flow rates at each point are established by heat transfer con- 
siderations alone. Several rigorous analyses for computing film thickness have been 
published (Ref 6, 7, 8 and 9). Test data have substantiated these equations. 

Accurate computation of the film thickness and interfacial energy dissipation depends on 
the designer's knowledge of interfacial shear and its dependency on the structure of the 
liquid-vapor interface which has as yet not been either experimentally or theoretically 
established. 

After the annular film thickness has been established, the stability analysis can be initiated. 
This analysis will determine whether the annular flow pattern indeed exists and if the pro- 
posed design will meet the objectives of pressure drop. 

No suitable general stability equation valid for all situations is available. The designer 
can either use flow regime maps - such as the map illustrated in Figure 15 - Page 33 
on which the local condensing flow conditions can be superimposed or he can locally apply 
adiabatic stability criteria involving the Weber, Reynold's, Froude numbers and the property 
ratios along the condensing film. A procedure such as this requires that the growth rate 
of the instability be large compared to the velocity of the liquid film. The effect on 
stability of the vapor velocity normal to the liquid interface caused by condensation is, of 
course, ignored in this approach. 

When instabilities manifest themselves, the designer should answer the following logical 
questions: 

1. Is the instability objectionable? 

2. Are objectionable low frequency, high amplitude liquid slugs formed, or is 

spray formed which may be tolerable? 

3. How is the heat transfer and pressure drop affected by the altered flow regime? 

This requires some knowledge of what is objectionable and what is tolerable. In addition, 
the designer must postulate the size of the detached and entrained liquid. This estimate 
could be scaled to wavelength. Unstable waves which are small compared to the tube diameter 
will form spray and those which are large compared to the tube diameter or vapor space will 
form slugs. The size of the liquid slug or drop and velocity of flow determine the frequency 
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and amplitude of the pressure pulses* 

In any event, the stability of wetting films must be predictable. The present state of the 
art is reviewed below and criteria for film stability are presented. 

4. 1.1.1 Review of Wetting Film Stability 

Most of the existing research concerning the breakup of liquid films was motivated by the 
practical problem of liquid film cooling for rocket nozzles. It is basically experimental, 
although some simplified analyses were also made. Much of this work was done for plane 
flows. This will not be extensively discussed herein because the basic configuration of 
interest for condensers is annular flow. Emphasis will, therefore, be given to the studies 
dealing with that type of flow. 

The striking aspect resulting from a survey of the literature on annular flow is that many 
different criteria are suggested for the film breakup and some work indicates that parameters 
found by others to be significant are not important. Relatively few attempts are made to 
explain the differences; the nature of most of the papers is the description of the particular 
studies made. For example, one of the most extensive experimental investigations was 
performed by Dunkler (Ref 10 and 11) who found that over the range of conditions which he 
studied in the downward annular flow in a vertical tube, there were essentially two signifi- 
cant changes in the nature of the flow. The first of these he associated with a basic change 
in the wave structure which occurs when the energy transmitted across the interface reaches 
some critical value and the second he thought to be associated with a liquid film velocity 
effect which decreases the wave amplitude. Dunkler further concluded from his results that 
the liquid Reynold’s number was not a significant parameter which determines wave motion 
and wave height. 

The experiments described by Kinney et al in Reference 12 were performed in horizontal 
ducts. They found that the liquid film surface became wavy when a critical liquid flow 
rate was exceeded. Furthermore, for more viscous liquid films this critical flow rate was higher. 
One can thus infer that in contradistinction to the work by Dunkler on a single-phase liquid 
film on the outside of a vertical cylinder, the liquid Reynold’s number is indeed a signifi- 
cant parameter. Brauer (Ref 13) was also led to believe from his experiments that the 
occurrence of surface waves is associated with the transition of the liquid film from laminar 
to turbulent flow, i.e., is related to a critical liquid Reynold's number. The experiments 
of Reference 12 covered a wider range of liquid viscosities than did Dunkler's experiments; 
this is perhaps an explanation of this discrepancy. It is also important to note that because 
the authors of Reference 12 seemed to be unaware of the basic governing parameters and 
physical mechanisms, their tests and correlations were unduly extensive, i.e., they varied 
the liquid velocity by itself, then the viscosity by itself, etc. 

More seemingly contradictory results were obtained by Knuth (Ref 14) who found surface 
waves for aJJI liquid flow rates in horizontal ducts. For liquid flow rates larger than some 
critical value he found a second type of surface wave (ones with long wavelengths). Thus, 
there now appear to be two distinct and different interfacial instabilities. The disagreement 
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between the results of References 12 and 14 was attributed to the differences in the methods 
of liquid injection in the experiments. Knuth, however, agreed with Kinney et al that the 
instability depended on the liquid Reynold's number. Recall this disagrees with Ounkler's 
conclusion. Further, References 12 and 14 are in agreement that the gas Reynold's number 
has a negligible effect on the interfacial stability whereas Laird (Ref 15) shows that over 
isolated ranges of operation the gas Reynold's number does influence the interface. Other 
work could also be cited to show the types of contradictions prevalent in this regard. 

Although the existing work seemingly contains contradictions and disagreements, it does 
indicate that objectionable phenomena can occur in two-phase flows when the gas-liquid 
interface becomes wavy. It is, therefore, essential to discuss the possible causes of such 
interfacia! waviness from a fundamental viewpoint. A body of fundamental knowledge exists 
which does not appear to have been generally known heretofore to the people investigating 
interfacial stability in two-phase flows. 

There are essentially four different types of instabilities which singly or in combination can 
cause an interface to become wavy and can thereby lead to film breakup or slugging in the 
flow: (1) hydrodynamic instability, (2) Kelvin-Helmholtz instability, (3) Rayleigh-Taylor 
instability, and (4) Benard instability. 

Hydrodynamic instability is the case in which a fluid is undergoing transition from laminar 
to turbulent flow as a result of amplification by viscosity of infinitesimal disturbances in the 
fluid. These disturbances originate either inside or outside the film. (Actually, all in- 
stabilities to be discussed here are hydrodynamic, but special meaning is given to this phase 
herein for convenience.) Hydrodynamic stability of boundary layers in a homogeneous 
fluid has been extensively studied (see Reference 16 for a summary), and it has been con- 
clusively established that such waviness occurs under certain conditions. 

The above hydrodynamic stability analyses were recently extended in Reference 17 for 
nonhomogeneous fluid flow. In particular, consideration is given therein to the specific 
case of two contiguous, viscous, incompressible fluids in plane motion; one fluid is bounded 
by a solid wall below and the other fluid above. The second fluid is unbounded above. The 
fluid motion is steady and unidirectional, parallel to the interface, and the shear rate in each 
fluid is uniform. As in the case of the stability of homogeneous fluids, the mathematical 
analysis is based on small disturbance theory and leads to an eigenvalue problem in a system 
of two linear ordinary differential equations. In addition to the dimensionless wave number, 
disturbance phase velocity, and the (inner fluid) Reynold's number occurring in single-fluid 
studies, the viscosity and density ratios and the Froude and Weber numbers appear as im- 
portant parameters for nonhomogeneous fluids. The effect of gravity is determined by the 
Froude number which represents the ratio of inertia to gravitational forces. The influence 
of surface tension is determined by the Weber number which essentially denotes the ratio of 
inertia to surface forces. Single-looped neutral stability maps (see Figure 1), similiar to those 
for homogenous fluids, and disturbance amplifications rates are presented in Reference 17. 

It is essential to understand that the work of Reference 17 is an extension of the hydrodynamic 
instability analyses made for single fluids and, hence, gives an indication of the conditions 
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SINGLE-LOOPED STABILITY MAP OF FELDMAN 
NEUTRAL STABILITY CURVES: WAVE NUMBER VERSUS LIQUID REYNOLDS NUMBER 
GRAVITY AND SURFACE TENSION FORCES ARE NEGLECTED 
LIQUID-GAS VISCOSITY RATIO = 10 
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under which the interface can become wavy as a result of amplification of disturbances in the 
liquid film bounding a rigid plane. Considerable confusion appears in the literature as to 
the type of instability studied therein (see Reference 18 for example). In essence the work of 
Reference 17 concerns the hydrodynamic stability of a film on a surface with new constraints, 
due to another fluid, at its edge. Furthermore, it must be clear that these results pertain 
only to plane flow of two fluids under adiabatic conditions with no mass transfer. Therefore, 
utilization of these results for condensing annular films (as are obtained in tubes) should be 
made with caution. 

The Kelvin-Helmholtz instability is the second possible type to influence two-phase flows. 
This instability arises when the different layers of a stratified heterogeneous fluid are in 
relative horizontal motion and is due to the interaction between the fluid media at the 
interface. The mathematical analysis of this instability is described in Reference 19. 
Basically it is also a small perturbation approach, but since attention is focused on the in- 
terface, the viscosities of the fluids are neglected. The Reynold's number, therefore, plays 
no important role in this type of instability. Instead a parameter (u^- uj 

for fluids with discontinous velocities and densities appears as the criterion, where g is the 
acceleration of gravity, a nd are density ratios defined as o< x - fi/f, 4 • , and 

U] and U 2 are the velocities of the two fluids. For fluids with continous velocity and 
density distributions in the vertical direction, this parameter is called the Richardson number, 
J, and is written as J =■-(£ This re P resen t s the ratio of buoyancy and 

inertia forces. For stability, neglecting surface tension, the disturbance wave number must 
be less than some value of the dimensionless parameter written in either of the above ways. 

The striking aspect of the Kelvin-Helmholtz instability is then that it occurs no matter how 
small the velocity difference or shear of the two fluids. The instability arises by the 
crinkling of the interface by the shear that is present. This crinkling occurs even for the 
smallest differences in the velocities of the two fluids; it can occur even if the motions of 
both fluids are laminar. The source of the Kelvin-Helmholtz instability lies in the energy 
stored in the kinetic energy of relative motion of the different layers. The tendency to- 
ward mixing and instability will be greater, the greater the prevailing shear force. The 
only counteracting forces are inertia and surface tension. 

It is important to note from the stability criterion cited above that in a zero-gravity en- 
vironment the two-phase flow would be unstable to all wavelengths. In Reference 20, a 
similiar criterion is derived by force balances for the Kelvin-Helmholtz instability. The 
effect of gravity has been omitted in this work. However, growth rates for this type of 
instability are determined which should be of value for the present work. (See Figure 2). 

The third possible type of instability is the Rayleigh-Taylor instability. This is an instability 
of the interface between two fluids of different densities which are stratified or accelerated 
towards each other. This type of instability arises from the character of the equilibrium of 
heterogeneous fluids. The mathematical analyses of this instability are also presented in 
Reference 19. They are either of a normal mode or a variational type. Studies have been 
made of both viscous and inviscid fluids, including and neglecting surface tension effects. 
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GROWTH RATE VS WAVE LENGTH FOR VARIOUS 
WEBER NUMBERS (REF 20) 



FIGURE 2 

The stability in the Rayleigh-Taylor sense depends only on the relative orientation of the 
density gradient between the fluids and the acceleration field. For example, horizontal 
fluid layers in a gravitational field will be stable only if the lower fluid is heavier than 
the upper; this configuration is unstable (neglecting surface tension) for all wave numbers 
if the reverse is true* For some unstable situations modes of maximum instability exist, and 
the dependence of the disturbance growth rate on its wave number has been explictly de- 
termined for a number of special cases. 

The last type of instability capable of influencing condenser flows is associated with having 
a fluid configuration in which the density variation is such that heavier fluid is above lighter 
fluid in a gravitational field. This is clearly a top-heavy situation which is potentially un- 
stable. The instability creates a tendency for the fluid to redistribute itself to remedy the 
weakness in its arrangement. However, the natural tendency of the fluid will be inhibited 
by its own viscosity. In other words, it can be expected that an adverse density gradient 
which is maintained must exceed a certain value before the instability can manifest itself. 
Considerable study has been given to this type of instability where the density gradient in 
the fluid results from a temperature gradient. However, similiar results will be obtained if 
the density gradient results from concentration gradients. This type of instability is called 
the Benard instability. The Rayleigh-Taylor stability applies to two fluids with an interface, 
whereas the Benard stability applies to a homogeneous fluid with density gradients. In a 
horizontal annular condensing flow the upper part of the fluid will be subject to an ad- 
verse density gradient as a result of both thermal and concentration gradients; the lower 
part will be in a stable configuration. The flow on the other parts of the annulus will 
have the gravitational force normal to the density gradient and, hence, will be subject 
to conventional natural -convection phenomena. Therefore, extreme asymmetry of the 
flow can result from such a configuration. Note that this discussion of Benard instability 
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relates to a fluid at rest. If a gas flow is superposed on such a fluid configuration, longitudinal 
or transverse vortex rolls result. A description of these is given in Reference 21. 

The normal mode and variational analyses of the B6nard instability are presented in Reference 
19. it is shown that a critical value of the Rayleigh number must be exceeded before in- 
stability occurs. (The Rayleigh number is the product of the Prandtl and Grashof numbers 
where the latter represents the ratio of buoyancy to viscous forces.) 

It is important to note that a particular type of instability under given conditions is generally 
more unstable to one type of disturbanceor wave than another. (Rayleigh, in Reference 21, 
says, "Some kinds of disturbances produce their effect much more rapidly than others.") 

For example, it is well known in hydrodynamic stability theory that for incompressible flows, 
two-dimensional disturbances (transverse waves) are more destabilizing than three-dimensional 
ones (oblique waves). See Reference 16. However, for compressible flows the three- 
dimensional waves are more destabilizing. Furthermore, for flows with body forces, 
longitudinal rather than transverse waves are the most destabilizing. Therefore, some study 
should be made to see whether each type of instability of interest is associated with distinct 
wave patterns. In this way it might be possible to identify the various types of instabilities 
or to determine the most harmful type. In the course of the experiments on film breakup 
described in Reference 22, two-dimensional, three-dimensional, and roll-waves were all 
observed. 

Now that each possible typeof instability has been described, the distinctions among them 
must be made and their relations to the problem of condensing flows must be established. 

First of all, note that the hydrodynamic and Kelvin-Helmholtz instabilities are fundamentally 
associated with fluids which are in motion; the Rayleigh-Taylor and B&nard instabilities can 
occur in fluids at rest. The primary factor in these is the body force action. The latter 
two types would not, therefore, be of consequence in a true zero-gravity environment. The 
hydrodynamic and Benard instabilities can occur in homogeneous fluids as well as in heter- 
ogeneous ones, whereas the Kelvin-Helmholtz and Rayleigh-Taylor instabilities are associated 
only with nonhomogeneous fluids. 

In a normal or reduced gravitational environment, all four types of instability are clearly 
possible. The Rayleigh-Taylor and B6nard instabilities most likely would not be so important 
as the other two types for films that have moderate motions because the body forces would 
be small relative to the inertia forces. The Kelvin-Helmholtz instability would most likely 
occur before the hydrodynamic instability because the former has been found to occur even 
if both fluids are in laminar motion. Near the back end of condenser tubes (and the front 
end of boiler tubes), the liquid layers meet and the flow velocity is very small. Under these 
conditions the Rayleigh-Taylor instability could be important. Work on this type of in- 
stability as well as studies of water-air, mercury-air interfaces at zero-gravity conditions 
(in a KC 135 aircraft) have been done at TRW and are reported in Reference 23. 

In any event, it would seem that determination of the instability most probable to influence 
a given flow first could be qualitatively obtained by comparing the various existing stability 
criteria. However, it must beremembered that all the above-described analyses were for 
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plane flows and none Included the effects of heat and mass transfer which are an inherent part 
of the condensing process. Therefore, at best, predictions made on this basis for two-phase 
annular flows are tenuous. It is, however, interesting to note that even despite the limita- 
tions of the existing stability analyses, they can be used to make some of the seemingly 
anomalous results obtained from the various experiments appear to be reasonable. For 
example, recall that in References 10 and 11 the two changes in wave structure were found. 
The first change, caused by the energy transmitted across the interface, certainly seems to 
be the Kelvin-Helmholtz instability; the second change, associated with the liquid flow 
rate, appears to be the hydrodynamic instability even though Dunkler did not think that 
the liquid Reynold's number was significant. In Reference 14 two different types of surface 
waves were also found. The first appeared for all liquid flow rates in consonance with the 
results of the Kelvin-Helmholtz stability analysis, and the second waves which occurred 
beyond a certain value of the liquid flow rate or Reynold's number were caused by the 
hydrodynamic instability. It is possible that Dunkler, who did not obtain the first type 
of waves until a certain amount of energy had been transmitted across the interface, had 
less disturbances in his experiments than did Knuth. Therefore, Knuth obtained the first 
waves for all liquid flow rates. These experiments tend to support the conjecture made pre- 
viously herein that the Kelvin-Helmholtz and hydrodynamic instabilities are the more im- 
portant pair for flows with at least moderate liquid velocities and that the Kelvin-Helmholtz 
instability appears first. 

The waviness found in References 12 and 13 after a critical liquid Reynold's number was 
exceeded must have been caused by the hydrodynamic instability. The reasons why no 
other wave structure corresponding to the Kelvin-Helmholtz instability was reported in those 
papers are not clear. Either the investigators were not looking for them or the disturbances 
level was below that required to cause such instability. Note in Reference 13 the flow was 
that of a single-phase liquid external to a duct and, therefore, less disturbances were im- 
posed on the liquid film. More careful analysis of the data of those papers might sub- 
stantiate that disturbance level was low. 

The conclusions of References 12 and 14 that the gas Reynold's number is unimportant as a 
prime interfacial stability parameter certainly are substantiated by the stability analyses, 
because this parameter does not appear at all. On this basis the results of Laird (Ref 15) 
which showed an effect in isolated regions of the gas Reynold's number remain questionable. 
It may be that the gas-flow turbulence level was responsible for this result. Further dis- 
cussion of this point will be made subsequently. 

In a zero-gravity environment the Rayleigh-Taylor and B4nard instabilities would not exist 
because there is no body force. It appears that the Kelvin-Helmholtz instability would be 
the dominant one, at least for film speeds below a rather large value (i.e. , below a large 
Reynold's number). The reason for this conjecture, again based on the existing analyses, is 
that this instability can occur even if both fluids are in laminar motion, whereas the hydro- 
dynamic type occurs when destruction of liquid laminar motion begins. Furthermore, the 
stability criterion for the Kelvin-Helmholtz type indicates that the flow is unstable for all 
wavelengths when g goes to zero. Whether, in fact, the situation is as serious as this cannot 
be determined without more careful study of this type of instability, including such effects 
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as heat, mass transfer, and viscosity. 


It is known from hydrodynamic stability studies for homogeneous fluids that cooling of the 
layer adjacent to a rigid surface tends to stabilize it. Therefore, in a condensing flow 
it would seem that the layer would be less subject to hydrodynamic instability. Whether 
this is also true for the Kelvin-Helmholtz instability remains to be investigated. The in- 
fluence of viscosity on the latter type, which has heretofore been neglected, might also 
be found to be stabilizing because it moderates the velocity difference at the interface. 

There is one additional word of caution necessary to be added with regard to the use of 
existing stability analyses and that is that they all refer to essentially unbounded flow 
in the direction normal to the surface. In effect, the flows treated are external ones 
whereas for the problem at hand the flows are internal. Although, as has already been 
pointed out, this "geometrical" difference may be important by itself, there is another 
aspect which may be of even greater consequence. In an internal flow one part of the 
liquid layer can produce disturbances to another part. Also in the main (gas) flow there 
exist such disturbances in internal flows (to a greater degree than in external flows) as 
turbulence, regular and random sound, temperature spottiness and the like. These dis- 
turbances may cumulate in different and nonlinear ways and the theory takes no account 
of these. Therefore, there is some question as to whether stability analyses as described 
above have any meaning for internal flows. The lack of correlation of hydrodynamic 
stability theory with transition data for flow in shock is clearly shown in Reference 24. 

Laminar internal flows are shown to become unstable at Reynold's numbers at least an order 
of magnitude lower than that predicted by stability theory. This is the opposite of what is 
meant by flow instability. The question of the validity of hydrodynamic stability analyses 
for normal internal flows is discussed more fully in Reference 24, 

The Kelvin-Helmholtz instability may also be affected by the larger disturbances associated 
with internal flows. 

4. 1.1.2 Prediction of Condensing Film "Breakup " 

The method for predicting the "breakup" time or distance from the condenser inlet to the 
point where liquid leaves the interface as a dispersed phase or where the liquid may bridge 
the tube as a plug is formulated below for an adiabatic liquid film model. 

Assume a liquid-gas interface to be perturbed by disturbances of all wavelengths at a multi- 
plicity of sources, and that the disturbance having the maximum positive growth ratedominates 
the interface. Since the maximum growth rate disturbance can originate at any point along 
the interface, it will be further assumed that only those formed at the furthest upstream 
distance will dominate. This assumption precludes the existence of standing waves which 
is reasonable since none have been observed in previous experiments. This is illustrated 
in Figure 3. 
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WAVE GROWTH ON LIQUID FILM 



Figure 3 shows a continous formation of maximum growth rate waves at L * O • LT 
represents the distance at which liquid escapes from the interface. Waves of lower growth 
rate or waves formed at sources where L> Q a re eliminated from the picture. This 
assumption is valid if the maximum growth rate has a "steep" maximum as shown in Figure 4, 
curve (a). If the curve is flat as in curve (b), then the assumption becomes dubious. 

GROWTH RATES IN LIQUID FILMS 


GROWTH 

RATE 



FIGURE 4 


A single wave formed at L s O 


will appear chronologically as in Figure 5. 
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PROPAGATION OF A SINGLE WAVE IN A STATIONARY FLUID 



Figure 6 shows a wave as it may propagate in reality. It is simpler mathematically (for 
linear problems such as that of stability) to make a Fourier analysis of such a wave and consider 
the individual Fourier components. It is, therefore, generally assumed in stability analyses 
that all the waves are composed of simple harmonic waves of the form: 

let ( L — C ©) 

v<L,Y,e)=f(Y)e (Eq „ 

As a result of this assumption all derived quantities, such as the velocities and amplitude, 
are harmonic and grow or decay exponentially. Equation 1 also defines the complex wave 
propagation speed C 5 +\C\ whose real part, Cr , is the physical velocity of pro- 

pagation of this simple wave. The imaginary part, C* is related directly to the growth 
rate as will be seen later. The wave speed, C , is measured with respect to the un- 
disturbed fluid; therefore, the wave velocity relative to a point fixed in space will vary 
with the velocity of the interface flow as shown in Figure 6. 

PROPAGATION OF A SINGLE WAVE IN A FLOWING FLUID 


(C R + U) 


It is clear from Figure 6 that: 

(Eq 2) 

If Equation 2 is applied to a forced convection condensing process in which U. x and Cr 
may be functions of length, it is more convenient to write Equation 2 in differential form 
as: 



For adiabatic fully developed flows the wave propagation velocity and interface velocity 
will not vary in the downstream direction, so the second term in Equation 3 will vanish. 
For simplicity, this case will be considered first. Therefore, Equation 3 becomes: 


L. = e(c R +a 2 ) 


dL = u 2 (^ R + l)dG = C R (i-f 0 (Eq 4) 

Since the ratio Cr/u.r appears as a parameter in Equation 4, it is necessary to obtain 
an estimate of its order of magnitude so that the "history" of a wave with a maximum growth 
rate (as illustrated in Figure 6), originating at l_— O and becoming harmful due to breakup 
(liquid entrainment) or flow plugging, can be determined. From Feldman's analysis of the 
hydrodynamic instability of a plane liquid film the ratio Cr/\^ z is plotted as a function 
of the property ratios of the two fluids in Figure 7. Figure 7 indicates that this ratio is 
approximately equal to 0.1. Therefore, Equation 4 can be written as: 


d L - u z d9 


(Eq 5A) 


If for other instabilities or flow conditions or configurations 
limiting form of Equation 4 is; 


IX 2 , /c^< 1 then the °ther 


dL=c R de 


(Eq SB) 


The maximum value of this ratio will have to be determined from analysis or experience for 
the situation corresponding to annular two-phase condensing flows. In all situations other 
than the limiting cases (given by Equation 5A and 5B), Equation 4 must be used. 

It is important to note that the time, Q , in Equations 4 and 5 depends on the growth rate 
of the wave, which is related to Cj . This relationship is explicitly determined as 
follows. 
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MAXIMUM WAVE CELERITY AS A FUNCTION OF PROPERTY RATIOS 



FIGURE 7 


t 


The amplification of the wave is described by the vertical velocity component of a point on 
the interface (Y = b) moving with the horizontal propagation velocity of the disturbance 
CL S For this case it can be found from Equation 1 that 


IZ-\« 


(Eq6) 


Differentiating Equation 6 yields 

If <5 denotes the liquid film height, then -4-& where B is the wave amplitude. 

Equation 7 can be written in terms of the wave amplitude as 
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B 



(Eq8) 


upon assuming fhaf the liquid film height, & , does not vary with time. Integrating 

Equation 8 and setting the integration constant equal to zero without loss of generality gives 


B 

B 


- C; 


or — C<c: dO 

B 1 


(Eq9) 


It is evident from Equation 9 that o( C| Is the parameter related to wave growth and, there- 
fore, to the stability of the film. 


Combination of Equations 4 and 9 yields 
dB _ **^1 J ^ 

e _ UzfVu.t 0 <£q io) 

If the limiting forms of Equation 5 had been used, the denominator term in parentheses in 
Equation 10 would not appear; in one case ( U.x/£ R << I ) C ft would appear in place of 
. In Reference 17, expressions are formulated for O (positive growth rates), 

C \- O (neutral stability), / and o< for maximum growth rate waves. These 

results for hydrodynamic instability are presented in Figures 7 through 9 for the case where 
gravitational and surface tension effects are neglected. (For the Kelvin-Helmholtz in- 
stability the results of Reference 20 could be used).* These figures together with Equation 
10 or its limiting forms are used to determine the breakup or plugging flow length. For 
other instabilities or configurations, of course, corresponding figures would have to 
be determined. However, for illustrative purposes the results of Reference 17 will be used 
herein with cognizance taken that these results apply explicitly to thin adiabatic liquid 
films with linear (Couette) velocity profiles. 


MAXIMUM GROWTH RATE FACTOR AS A FUNCTION OF PROPERTY RATIOS 


*See Figure 2 



FIGURE 8 
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For adiabatic films the abscissa of Figure 8 is 


r u, r 


w*a> 


where 




4 

Integration of Equation 1 1 yields 

C-, = u 2 F^„ Re 4 + K 


where K is an integration constant. When Cj = 0^ Re* s R^ where 
the Reynold's number for neutral stability. Therefore n 

K - Re* 


and 




If Equation 12 is substituted into Equation 10 there is obtained 

A3 (Ref - RerJ) J 1 

B (Vu, + 1 ) 


Replace oc * by its dimensionless form used in Figure 9, 0<. - 

Equation 13 becomes 


45 = £L F w ^Cr^-R^ Al 

B ( Cr /ol 2 t I) ^ 


(Eq 11) 

Re 4« Is 

(Eq 12) 

(Eq 13) 

2H^ so that 

7\ 

(Eq 14) 
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MAXIMUM GROWTH RATE WAVE LENGTH AS A 
FUNCTION OF PROPERTY RATIOS 



FIGURE 9 


To find the breakup or plug flow length. Equation 14 must be integrated. If denotes 

the point where the waves lead to unfavorable phenomena such as^liquid entrainment and 
flow plugging, then the integration limits are L_ = 0 to L= L. .^Similarly, 3 s B® 
to 8 s 13* where is the initial wave amplitude and B is the^jmplitude at 

which either entrainment or plugging occurs. To determine explicitly 3 , the 

amplitude fortearing of the interfacial waves can be obtained from experiments, whereas 
for plugging it is clear that the wave amplitude plus the film thickness must be of the order 
of one-half the flow passage diameter. 

For the case of liquid entrainment it has been shown in References 22 to 25 that the wave 
heights approach the magnitude of the average liquid film thickness. Photometric measurements 
of the wave profile reported in Reference 25 are reproduced in Figure 10. This figure indicates 
that near the point of liquid entrainment B/^ = 0.81. Another such figure derived in 

a similiar way (Figure 4c of Reference 25 shows that B / £ = 0.974. 


Figure 1 1 taken from Reference 25 shows maximum net wave height measurements (2B) 
versus film thickness. These measurements are limited to conditions for which no liquid 
is entrained. The upper part of the curve represents the point of incipient entrainment at 
which 
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WATER-AIR WAVE PROFILES (REF 26) 
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FIGURE 10 


AIR-WATER MAXIMUM WAVE HEIGHT VERSUS FILM THICKNESS (REF 25) 



FIGURE 11 
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Other evidence of this sort is also available. It should be noted that these data were taken 
for the annular downflow in a vertical tube, and, hence, are more appropriate to the present 
work. However, these data were obtained for the liquid film application to rocket nozzles; 
therefore, the range of parametric values may not be analogous. 

For the other end of the integration range it is necessary to know the initial wave amplitude 
B q . Reference 27 states that ly\ &Vb 0 reaches a constant value for the breakup of 
moving liquid sheets. (This configuration is again not the same as the one of interest 
herein, but is used merely to obtain some orders of magnitudes from existing data.) 

Reference 28 theoretically predicts values for this ratio from which "B 0 can be determined 
when £&** is known. 

For the other undesirable case of flow plugging, assume that flow plugging will occur when 
“ & . Therefore, for this case the upper integration limit is 



which will occur when L - L s . If this length is less than that for entrainment, 
ijj* , slugging will actually take place in the flow passage. Alternate criteria for 
slugging may also be developed. For example, if the size of the entrained droplets as 
computed from Reference 29 are of the order of passage dimensions, then slugging can also 
occur. 

Now that the integration limits have been established. Equation 14 can be integrated to 
give 


\y\ §* - _ fcwjocCRe-f L 

( < ~ r /vjl 2 4* I ^ ~cS 


(Eq 15) 


Since the concern here is with the maximum growth rate wave, the maximum values of 

ar *d CR/Uashould be inserted in Equation 15. These can be evaluated from Figures 7, 
8 and 9 which, of course, are limited to a specific configuration and instability type as 
discussed previously. Accord ing to these figures Fm^x^max ant ^ /C^^^are 

functions of property ratios alone; therefore, for a given combination of fluids at a certain 
pressure and temperature these quantities are constants. Also, as was previously indicated, 
the ratio B / 0^ reaches a constant value for liquid film breakup. 


Equation 15 can be rewritten as 

L_ max G 




(Eq 16) 
(Eq 17) 
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where El is a constant for a given set of property values and 

and is explicitly defined as 


E> 


- »*> b */b. I( Cr A^ vry&x~^~ 0 


c* 




vm?yx 


(Eq 18) 


A 

To obtain an idea of how Equation 17 can be used to find the breakup length L and to 
indicate the data required to determine this length for the problem under consideration, use 
will be made of all existing data even though the configurations and ranges of conditions 
of the various sources are not consistent. Thus, if air-water systems are considered at 
atmospheric pressure and temperatures with ~ 54. 1 and ^y/^: 5=5 0.0012, 

then from Figure 7, for hydrodynamic instability of thin plane film, 
and = 0.6 from Figures 7 and 9. Note that this corresponds to a maximum 

growth rate wavelength-to-film thickness ratio of about 10 which corresponds to the visual 
experimental observations made in a tube with annular flow reported in Reference 30. From 
Reference 28 |\^ ^ and from Figure 8, F 

magnitude of E is 


_o 

= 4 x 10 . Therefore, the 


F - + Q 

11 ~ o.<b ( 4 x IO“*) 


- 5.5 * io 


This calculation is for film breakup for which liquid drops are entrained. As already pointed 
out, the integration limits (and therefore B /j^ 0 ) would be different for slugging 
conditions. As a check on the above calculation note that Equation 17 indicates that 
L*/£ is a function of . Data are plotted in terms of these quantities from 

eight different sources (primarily for annular flows) in Figure 12. Table I describes the 
test conditions in detail. Figure 12 indicates that the data lie approximately on a 45 degree 
line. The film thickness, e> , therefore, plays no role in the breakup length. This is 
consistent with the model used in deriving Equation 17. Thus, Equation 17 can be empirically 
adjusted to the form = 1.6 x 10^. If it is assumed that Rg>^, then 

Equation 17 can be written as 


'- xj ^ £ ± - - E = S.S x IO 9 

The value (5.5 x 10^) predicted by the method developed herein does not agree with that 
indicated by experiments (1.6 x 10^). The lack of agreement between the predicted and 
experimental results should not be too disturbing because the predictions were calculated 
essentially from the work of Reference 17 for hydrodynamic instability of plane adiabatic 
thin films in Couette flow. In the experiments the flows were mainly annular with heat 
transfer in some cases, and the Kelvin-Helmholtz instability could also have influenced the 
film breakup. Further, the analysis of Reference 17 which was based on the theory of small 
perturbations loses accuracy as wave amplitudes become sufficiently large. If the disagreement 
between predictions and reality persists after more appropriate data for both calculations 
and experiments have been obtained, the predictions would have to be based on a higher- 
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BREAKUP LENGTH VERSUS REYNOLDS NUMBER 



FIGURE 12 
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FILM BREAKUP TEST DATA 
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order theory to account for the fact that the amplitudes near breakup are large. Some study 
of the growth of waves with large amplitudes has already been made (Reference 31) and It 
has been found therein that when the wave amplitude has grown sufficiently It remains constant, 
i.e., the growth rate goes to zero. If this w^je so in the present problem. It would imply 

L (experimental) for fromj 


that 


L* 


(predicted) should be less than 


Equation 17 it can 


be seen that if (which is related to the growth rate) decreases, L"' increases. 

Clearly this is another aspect of the problem which should at least be kept in mind. 


The calculation method described herein appears to be more plausible, however, if com- 
parison is made with results more closely analogous to the conditions of the theoretical 
model. For example, for the special case of a single-phase boundary layer =1 and 

. Then according to Reference 32 for laminar-turbulent transition on a flat plate 




- 2 ..% * IO 
'+v-&Yisitio\n 


This number is, of course, of the same order of magnitude as that experimentally determined 
herein, as is to be expected. This also further substantiates the previous discussions that 
indicate the work of Reference 17 is directly related to hydrodynamic instability. 


In summary, it can be stated that the above derivation is general, i.e., it can be applied 
to any configuration and instability type. It was applied to a specific problem herein merely 
for illustrative purposes. To extend the applicability to non-adiabatic flows such as are of 
interest in condensing problems, the second term of Equation 3 must be dealt with. However, 
on the assumption of fully developed flow, this term can be shown to be of little significance 
and, again, only the basic data such as Cp/^^etc. must be found for the application 
of the method to this case. 


By applying the growth rates of Reference 20 (see Figure 2) to a mathematical treatment 
similar to the one already presented, it can be shown that the Kelvin-Helmholtz instability 
can be expressed in terms of a Weber number based on breakup length, namely, 

(T 

If we are consistent with our original premise that the maximum growth rate wave dominates 
the interfacial breakup, then we canassume that whichever instability, hydrodynamic or 
Kelvin-Helmholtz, gives the smaller breakup length, that will be the one of interest. This 
is rational since the greater the growth rate, the smaller the length. These assumptions are, 
of course, subject to experimental verification. 

4. 1.1. 3 Stability of Condensing Films 

Condensation introduces additional complexities into the mechanics of film stability. Con- 
densation requires some modification in the boundary conditions used in the development of 
the stability equations of References 17 and 20. Also & , VJ ^ and U 2 are 
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BREAKUP REYNOLDS NUMBERS VERSUS CONDENSATION 
VELOCITY RATIO 



functions of length L. for the case of condensation and are thus considered variables 
during the integration for breakup length. 

The boundary equation for the vapor-liquid interface contains an equality of normal, shear, 
and surface tension forces. Condensation contributes an additional shear and normal stress 
term created by the vapor velocity normal to the interface, , which depends on the 

magnitude of the heat flux. These stresses are, 


Condensation normal stress: 

■ , * *• 

AP„« 


(Eq 1?) 


(Eq 20) 


Condensation shear stress: 

T m c A U** 'Ay 

To determine the effect of these quantities on the growth rate factors ( C* or S ) will re- 
quire detailed analytical investigations. 


The effects of condensation on film stability can be investigated empirically by plotting the 
"breakup" Reynold's number versus the ratio UL /tjL v aS ' n Fig ure 13* When this quantity 
approaches zero, we are near adiabatic conditions. 


*/A “ ^A h v * 
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A designer desiring to apply the "breakup" Reynold's number to a condensing film must 
determine the interface velocity, U.^ . Normally the vapor velocity, U»v , is the 
known quantity. The ratio /lAv / the mean film velocity to the mean vapor velocity, 
is called the slip and is determined by an involved momentum transfer analysis which will 
not be covered in this report. The results of two investigations will be cited which should 
suffice, namely: 


Reference 21: 

*£ 

U. v 

Reference 24: 



(Eq 21) 


Uv “ \% Zck'" 

HI 1 

where o< is the volume fraction occupied by the 
two-phase flow systems under consideration. For thin 
assumed. Thus, U.^ = U. 2 / 2 . 



(Eq 22) 


vapor which approaches one for most 
films a linear velocity profile can be 


4.1.2 Horizontal Two-Phase Flow 

Two-phase phenomena are generally complex and this is more the case when heat transfer 
is also considered. It is possible, however, to make certain general and specific observations 
for the case at hand, namely, the forced convection wetting condensation of mercury in 
horizontal tubes. For this case, the general adiabatic two-phase flow observations and data 
should be applicable. In particular, flow regime can be identified and the appropriate 
pressure drop and inventory correlations applied. Such anapproach was successful for the case 
of Refrig erent-22 condensing in horizontal tubes (Ref 33). In this instance, the Martinelli- 
Nelson approach was used to predict pressure drop to within ±30%, Similar ability to predict 
pressure drop is expected for most of our range of operations as the primary regime anticipated 
is that of annular flow. This regime should be amenable to the Martinell-Nelson equations. 

In terms of the fluid mechanics problem, the effect of a heat transfer is to cause an apparent 
vapor velocity normal to the wall. This has two consequences. There is a momentum change 
with length for the fluid passing through the pipe and, in addition, the flow configuration is 
altered as compared to the same flow rates for each phase without any heat addition. Because 
visual observations are so difficult when there is heat addition, there are very few such 
observations. 

One of the most striking characteristics of a two-phase flow is that, in general, the phases 
do not move at the same velocity. The gas phase usually moves more rapidly than the liquid 
phase. Because of this, it has been found convenient to speak of a slip velocity ratio. If 
this ratio is specified, then it is possible to determine the density of the mixture in the pipe 
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from the flow rates of each of the phases. The slip velocity ratio is related to the density 
in the following manner: The continuity equation on each phase yields 


tu v - lu t X - 


U.y Av 
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(Eq 23) 


A v and A f are the cross-sectional areas of a pipe occupied by the vapor and a liquid, 
respectively. If f \ p is the cross-sectional area of the pipe, then ^v/Ap » s fraction 
occupied by the gas. The static quality X 5 is the ratio of the weight of vapor in a section 
of pipe to total weight of liquid and vapor in that section. For an arbitrary length of pipe, 
the static quality becomes 
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When Equation 23 is substituted then 


v - iXJ ° 
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Thus, on re-arranging, 

X _ Uv / *S \ 

i -x ” u f \ i -x s / * 24 > 

Particularly at low pressure, the vapor phase normally moves with many times the velocity 
of the liquid phase, so the static quality is much lower than the flowing quality. 

Several different experimental methods exist for determining the static quality and there- 
fore the slip. One consists of suddenly isolating a section of pipe in which the flow is taking 
place and measuring how much liquid is present. Another method is to use a collimated 
gamma ray source and measure the attenuation in passing through the pipe. If the flow con- 
figuration is simple enough, a photograph can often yield the desired information. 

When two phases flow in a horizontal pipe they can be distributed in the pipe in a number 
of different configurations. These configurations are called flow regimes and the analysis 
of any two-phase flow problem begins by specifying what the flow regime is. 
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Figure 14 illustrates the flow regimes that can exist in a horizontal pipe while Figure 15 
shows a flow map that allows one to determine what the flow regime is when the pipe size 
and the flow rates of the two phases are know. The data of an additional reference (Ref 34) 
has been placed on this curve. It should be stressed that the divisions shown as lines in 
Figure 15 should really be bands. In addition, it should be stressed that this map applies 
only In fully developed flow for a wetting fluid without heat addition. The coordinates for 
this map are the ones suggested in Reference 35, while the data used to place the divisions 
are largely from References 36, 37 and 38.’ Insofar as most of the data were taken on pipes 
approximately one inch in diameter with air and water at low pressure, the placing of these 
lines must be regarded as tentative. Several of the flow regimes illustrated in Figure 14 
are combined in the map of Figure 15. This has been done because the same analysis for 
the pressure drop and density applies for each of the groups of flow regimes. 


HORIZONTAL PIPE TWO-PHASE FLOW REGIMES 



Stratified 



Froth Annular 



Mist 

* For critical flow FIGURE 14 


Fairly complete empirical analysis of flow in large horizontal pipes is reported in Reference 
37. 

The major portion of any condensing rig Is in the annular region but, depending on the tilt and 
the discharge end connections, it may go into the slug or stratified regime. To have stratified 
flow, of course, the tube exit quality must be greater than zero. A typical condensing run is 
plotted in Figure 15, which shows that annular flow is expected from 100% to about 5% 
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HORIZONTAL PIPE TWO-PHASE FLOW REGIME MAP FOR 
ADIABATIC WETTING FLUIDS 
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FIGURE 15 
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quality or until the vapor velocity is less than 10 ft/sec. (see Appendix C). For theconstant 
diameter tube, this will not occur until the last several Inches of the condensing section. 

At this point, stratified and/or slug flow would exist, depending on end conditions. Some 
pressure pulsations could occur even though the pressure drop in this region is expected to 
be essentially insignificant. 

End conditions might determine the flow regime if one has an exit quality greater than zero, 
since the condensing run passes between the two slug-stratified transition limits shown, one 
transition being for a filled exiting pipe (Ref 38 and 39) and the other for an unfilled exiting 
pipe (Ref 34). 

4. 1.2.1 Annular Mist and Foam Flow in a Horizontal Pipe 


Martinelli and others (Ref 40, 41, 42, 43) have developed the most general correlation 
available for calculating the pressure drop and density in a two-phase flowing system. Though 
there are some difficulties and internal inconsistencies in the analysis, it can be said for most 
of the region in which it is supposed to be valid that the Martinelli method predicts the pressure 
drops to within ± 30%. In the remainder of this section, the Martinelli method will be ex- 
plained and its limitations pointed out. 

Martinelli found in his original experiments that the pressure drop in two-phase flow is 
essentially independent of the details of the flow configuration so long as: 

1. the flow is steady, and 

2, there are no radial gradients in pressure in the tube, i.e. , the static pressure 
drop in the liquid is equal to the static pressure drop in the vapor. 

Essentially It is assumed, then, that an annular, annularsprayorfoam-flow pattern exists. Slug, 
wave and stratified flows are excluded from the pressure drop correlation. 

To begin, Martinelli separated the various possible two-phase flows into four mechanisms. 

These mechanisms do not correspond directly to any particular regimes of flow but depend 
on whether the Individual phases flowing in the pipe would be viscous or turbulent. They 
really should be regarded a$ a correlating device. These flow types or mechanisms are as 
follows: 

1. Flow of both the liquid and the gas may be turbulent (turbulent-turbulent 
flow). 

2. Flow of the liquid may be viscous and flow of the gas may be turbulent (viscous - 
turbulent flow). 

3. Flow of the liquid may be turbulent and flow of the gas viscous (turbulent-viscous 
flow). 
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4. Flow of both the liquid and the gas may be viscous (viscous-viscous flow). 


From the experimental data available at the time, the dividing line between the viscous 
and turbulent regimes was set as follows: 


Flow Mechanism 

Symbol 

Re fp 

Re 

— SE 

T urbu 1 ent-turbu 1 en t 

t-t 

> 2000 

> 2000 

Viscous-turbulent 

v-t 

< 1000 

>2000 

T urbu 1 ent-v iscous 

t-v 

>2000 

< 1000 

Viscous-viscous 

v-v 

<1000 

<1000 


The region 1000 K Re < 2000 is a transition region. To calculate the pressure drop 
which would take place in a two-phase flowing mixture, a quantity was defined which re- 
lated the two-phase pressure drop to the pressure drop of either component flowing alone at 
Its mass flow rate in the pipe. The definitions are as follows: 

4>*= f A % L ) T p/^AP/a l) ^ / 

§£- (^/ALW /(AP/al)* * 

The liquid and vapor fractions are deFined below: 

fraction of pipe area occupied by liquid, 

|^ = fraction of pipe area occupied by gas. 

The four quantities are just those needed to predict the pressure 

drop and density of a two-phase flowing system. By physical reasoning Martineili postulated, 
subject to experimental verification, that they are all functions of a single variable which 
is defined as follows* 



This ratio can be calculated from the mass flow rates of the two components and the fluid 
properties. When this is done the result is 




(Eq 26) 


The exponents vn 


and V\ and coefficients 


Cf and C. A are presented below: 
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v-t 
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v-v 

v\ 

0.2 

1.0 

0.2 

1.0 

Vrv 

0.2 

0.2 

1.0 

1.0 

C* 

0.046* 

16.0 

0.046* 

16.0 

c * 

0.046* 

0.046* 

16.0 

16.0 


These values for the exponents and the coefficients are determined from the appropriate 
single-phase expressions for laminar and turbulent flow. Figure 16 shows & and . 
flow mechanisms plotted on the same curve. o V 


Only one R$. curve, presented In Figure 16, was found to correlate the holdup data for 
all mechanisms. Later, in attempting to apply these results to steam-water mixtures at 
high pressure (Ref 42), it was found that the pressure drops were definitely over-predicted 
Because of this, was allowed to take on an empirical pressure dependence. The low 

pressure curves apply to $4 4-4 * I* 1 principle, they can be adapted for mercury 

condensing, but it is not recommended tnat they be used for the following reasons: 


1. A large part of the region of interest is viscous-turbulent; for this region there 
are no pressure-corrected curves. 

2. The two-phase pressure drop is tied to the single-phase liquid pressure drop. 

This view is more appropriate for boilers where 100% liquid enters and very 
low quality steam leaves than for condensers where vapor enters and most of the 
pressure drop occurs in the high-velocity range, vapor-fraction region. 


If one examines the pressure sensitive .. curves of Reference 42 clos 
seen that most of the pressure sensitivity could be eliminated by using ^ 


curves of Reference 42 closely, it can be 


In the light of these considerations, it is recommended that Figure 16 be used along with a 
numerical integration up the tube. 


* For smooth pipes 
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$ 9 versusXfor the viscous-turbulent and turbulent- 

turbulent MECHANISMS AND 'RfVERSUS'X (FROM REF 40) 



FIGURE 16 


TWO-PHASE SLUG FLOW FRICTION FACTOR MULTIPLIER TO BE USED IN EQUATION 27, 

(FROM REF 34) 
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FIGURE 17 





The original curves of Martinelli wer-e taken at low pressure and should apply quite well 
to the case of mercury in the temperature range of 600° to 750°F. The values of 
which are appropriate are very small, so it is recommended that the 5^ curves be used 
because the variation of JE with X in this region is very small. In fact, a constant value 
for is probably adequate, 

4. 1.2. 2 Slug, Plug and Wave Flow 

This flow regime is analyzed in Reference 44 and the method proposed is based on that 
reference. For a straight tube the pressure drop in this regime is probably negligible. 


Pressure drop is calculated from Figure 17 and the following equation: 
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ap tp = 


f-5 


(uJt + UJgf 

A. 


(v,+v^W 


•f = Moody friction factor 

multiplier from Figure 17. 


(Eq 27) 


The above equation is the usual pressure drop equation containing an arbitrarily defined 
mixture density. 


The Reynold's number for obtaining the friction factor is defined below. The friction factor 
f can be obtained from the usual friction factor curve using this Reynold's number 
definition: 



TT Dp 


4 ( tUf + 



» 


(Eq 28) 


To evaluate the density in slug flow. Reference 44 suggests that (the fraction of the 

pipe filled with gas) for an absolutely horizontal pipe is 



(Eq29) 


This expression works well over the entire range of conditions in the slug, plug, wave flow 
regime of the flow regime map. It fails only when the pressure is close to the critical where 
the slug flow hardly exists anyway. 
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4. 1.2. 3 Stratified Flow 


At very -low -volume flow rates the liquid and gas phases stratify with a smooth separation 
between them. As a rule, the pressure drops are very small . It has been found that usually 
the gas shear stress on the upper surface of the liquid is not very important in moving the 
liquid through the pipe, but that gravity plays the most important role. This means there is no 
"fully developed" flow but that for adiabatic flow the liquid level drops as the end of the 
pipe is approached. Thus the pressure drops are of the order of the tube diameter. The 
height of the liquid in the pipe and the pressure drop are both functions of the liquid of the 
pipe except for very long pipes. 


Stratified flow is of no significance for the case of zero gravity condensing. However, tests 
must be run on earth and understood since a significant amount of inventory may be involved. 


The most complete description of stratified flow is given by Bergelin and Gazley(Ref 38). 
Measurements of gas pressure drop and fluid heights were made in one-and two-inch diameter 
tubes, 10 and 16 ft long, for an air and water system. The gas -phase pressure drop data for 
the two-inch pipe was correlated by means of an apparent Reynold's number and friction 
factor. The friction factor is defined in the following equation! 


/4 p\ _ 4f ^1 

^ic/W D * 


while Reynold's number is defined as 

__ 4 £ 

” tt D P fig (Eq3,) 

A typical plot of friction factor versus Reynold's number is shown in Figure 18* Figure 19 
gives the liquid heights in the pipe ten feet from the. free overfall. 

The friction factor curves of Figure 18 are somewhat misleading as the true area for gas flow 
is not the pipe area. If the measured liquid heights are used to find the true area available 
for gas flow and the M hydraulic diameter" concept is used to compute a Reynold's number, 
the friction factor versus Reynold's number curves all fall very close to that for a smooth 
pipe in both the laminar and turbulent region. When a sufficiently high Reynold's number 
is attained in the liquid, the interface becomes wavy and the friction factor increases very 
rapidly. A large fraction of the total liquid in the system is probably in that part of the 
system which might be in stratified flow. The typical run shown in Figure 15 indicates that 
whether the nappe springs free at the end or not and whether the tube runs full depends on 
the end connections. The difference between the full and not -full case is quite substantial. 
It is also not clear whether liquid inertia and vapor shear will substantially assist in carrying 
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FRICTION FACTOR VERSUS REYNOLD'S NUMBER FOR STRATIFIED FLOW 
IN A 2.065 INCH ID PIPE WITH VARIOUS WATER FLOW RATES 


Lb/Hr 


The friction factor is defined 
by Equation 30 and the 
Reynold's number by Equation 
31 . These data are really In 
an entrance region. (REF 51) 


4-0 log (Pci? 


FIGURE 18 






Interfacial Height 10 Feet Upstream from Free Overfall, Inch 


INTERFACIAL HEIGHT IN A 2.065 INCH ID PIPE WITH CO-CURRENT FLOW 
OF AIR AND WATER (REF 51) 
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FIGURE 19 



PRESSURE FLUCTUATIONS IN HORIZONTAL SLUG FLOW (FROM REF 34) 
CROSS-PLOTTED FOR SMALL DIAMETERS 



FIGURE 20 
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REDUCED FREQUENCY OF PRESSURE PULSATIONS AS A FUNCTION OF THE GAS 
CONTENT FOR THE FLOW OF WATER AND AIR IN HORIZONTAL PIPES 
WITH DIAMETERS FROM ABOUT 1 TO 4 INCHES (REF 34) 



FIGURE 21 

SPEED OF MIXTURE UmAT THE STRATIFIED FLOW LIMIT AS A FUNCTION 
OF GAS CONTENT AND PIPE DIAMETER (FROM REF 34) 
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FIGURE 22 
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the liquid along in stratified flow. 

In any case, it is desirable to have a method for estimating the density in stratified flow. 

Figure 19 gives this. Once stratified flow occurs, the depth is quite constant as a "dam" at 
the discharge limits the flow. Chato (Ref 39) found a short ramp of liquid at the beginning 
of stratified flow leading up to the virtually constant depth portion. This ramp is largely 
supported by the shear of the vapor. 

The depth from Figure 19 is the depth at the discharge section. Equation 32 from Reference 
39 for an air-water experiment or a rough calculation going back from the discharge will 
give the depth inside the tube, it is not recommended undertaking any calculation here 
though, unless the experiment indicates stratified flow is an important flow regime. 

In Reference 39, flow regime observations were made for inlet velocities of 6 ft/sec or less 
and it was found that the conditions of discharge from the tube were very important in 
determining how much inventory there was in the tube. Tilt was also found to affect inventory 
tremendously, with one degree change having a significant effect. 

4. 1.2. 4 Pressure Fluctuation Frequency and Amplitude 

Both these quantities can be estimated from the adiabatic information given in Reference 34. 
Significant fluctuations occur only in the slug, plug, and wavy region. One will usually 
pass through this region, however, in the last few inches of the condenser, so the fluctuations 
may never develop any great amplitude. As a result, the results presented below will be 
conservative. In addition, the whole range of fluctuation amplitudes will (theoretically) 
be passed through so that one hardly knows what to expect in an experimental measurement 
in a condenser. Some of the data given in Reference 34 were located on the flow regime map. 
Figure 15, and were found to lie in the middle of the slug flow regime. The data of Reference 
34 were taken without heat transfer in 1, 2, 3 and 4 inch pipes with air and water at low 
pressures. Pressure oscillations amplitudes, and frequencies were measured and correlated in the 
the manner shown in Figures 20, 21 , and 22. 


Pressure fluctuation amplitude may be estimated from Figure 20 where 




is the mixture velocity = 


Yi±_ v t 

A P 


is the stratified slug transition velocity and a function of pipe diameter 
and quality as given in Figure 22. 


a is the fluctuation double amplitude for water. Because of the relationship 

of pressure drop to liquid density at constant velocity (AP*^ 3 ),it is re- 
commended multiplying this value of a. by a constant faijp to estimate 
the pressure fluctuations with mercury. ® (H^O 


Figure 20 was obtained by cross plotting and extrapolating the larger tube size data of 
Reference 34 to the smaller tube size range of interest. 
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The Frequency of pressure fluctuation can be obtained from Figure 21 taken from Reference 
34 directly. In Figure 21 , the actual frequency, , is obtained for a given gas content 

via the reduced pressure parameter N . and U. m are defined above. 

4.2 Forced Convection Condensation of Non-Wetting Mercury in Zero Gravity 

A drop-type flow regime is expected during the non-wetting condensation of mercury in 
zero gravity or with high vapor velocity. The condensed mercury is expected to condense 
on the surface and form a thin film (approximately 1 to 10 microns) which grows until it 
becomes unstable and breaks into drops (Ref 45). The drops formed have been observed in 
transparent tubes (Ref 4 and 5) and via the fluoroscopic observations of this program. The 
drops were observed to grow to a critical size and be accelerated toward the interface by 
vapor flowing over them. 

The non-wetting property and high surface tension often encountered with mercury rarely 
result In stable films of pure liquid mercury in a flowing system. Droplets dispersed In a 
vapor phase was the flow pattern observed In all two-phase flow tests made In transparent 
tubes using non-wetting liquid mercury In either nitrogen gas or in mercury vapor (Ref 5). 

The non-wetting flow regime map of Figure 23 resulted from these tests. 

In postulating either a heat transfer model or a fluid dynamic model for analytic purposes, a 
flow pattern involving liquid drops is therefore assumed. In condensing it is assumed that 
growing stationary surface drops are swept into the flowing vapor when they reach a certain 
critical size. The critical droplet size is determined by the vapor velocity, as shown in 
Reference 5. Since gas velocity will usually vary axially in a condenser, droplets of 
different sizes result. Figure 24 illustrates the flow pattern encountered during condenser tests. 

The basic flow regime for the non-wetting'condensation of mercury is then one of droplets 
on the wall and in the vapor stream. Droplet distribution is, of course, affected by gravity. 

In zero gravity a pure random drop distribution created by the mechanism described above 
is expected. With vertical upward and downward condensation, a significant effect on the 
phenomena would be expected. For example, in vertical upward flow higher velocities are 
required to prevent drops from falling opposite to the direction of flow. Vertical condenser 
design considerations were reviewed by TRW in Reference 5. 

The two-phase static pressure gradient of a condenser with a flow pattern as depicted in 
Figure 24 can readily be obtained from the momentum equation. The two-phase static 
pressure gradient consists of three terms, namely. 
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! (Eq 32) 
momentum 


The most difficult term to evaluate is the first term since it must encompass drops of all 
sizes. This can be accomplished by means of a distribution function as follows: 
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FLOW REGIMES FOR NON-WETTING MERCURY 
FLOWING WITH NITROGEN IN GLASS TUBES 
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Drop Flow Turbulent Flow 
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FIGURE 23 


NON-WETTING MERCURY CONDENSING FLOW PATTERN SHOWING: (1) DROP ON HEAT 
TRANSFER SURFACE GROWING DUE TO CONDENSATION, (2) DROPS IN TRANSIT FLOWING 
TOWARD SUBCOOLER INTERFACE, (3) DROPS OF DIFFERENT DIAMETERS DUE TO A VARIABLE 

AXIAL VAPOR VELOCITY 
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Typical Drop Trajectories 


Subcooler Interface 


FIGURE 24 


46 




dp 

dL 


drag of 
transit drops 



£ C (U 


(Eq 33) 


where N is the distribution function and N represents the number of drops per unit 

volume having diameters between and + <!& at distance L_ along the tube. Equation 
33 is derived in Appendix B, where N is shown to be 


M = ' 
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(Eq 34) 


The term (- is evaluated by tests as described in Reference 5. For zero gravity and 

horizontal tubes , 



(Eq 35) 


(Eq 36) 


where E 0 - is a constant as evaluated by experiment and represents the amount of defor- 
mation from the action of surface forces (drag) and body forces that a stationary drop can 
withstand before it is entrained and transported by the flowing vapor. 

The quantity € in Equation 33 represents the velocity ratio of the liquid (€ = and 

is determined by means of the equations of motion as covered in Appendix B. In order to 
evaluate £ the following non-linear differential equation must be solved*. 


du - 4 7 ^ “3T € X dL (Eq 37) 

(l-€)| l-€| is used rather than to account for proper sign in case 

which will occur in a constant diameter condenser. The coefficient of drag £4 can be 
evaluated by means of Figure 25 which takes into account the deformation of liquid drops 
due to drag. 


The second term of Equation 32 is the pressure gradient induced by the "rough" surface 
composed of the wall with stationary growing (condensing) drops. An equation of the form. 
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(Eq 38) 
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DRAG COEFFICIENTS FOR SPHERES, DISKS AND DROPS * 



FIGURE 25 


FRICTION FACTOR FOR SMOOTH AND SAND ROUGHENED 
PIPES AFTER NIKURADSE (REF. 48, FIGURE 84) * 



FIGURE 26 


*See Hughes, R.R. and E.R. Gilliland, "The Mechanics of Drops" 
Chem. Engineering Progress . Vol. 48, No. 10, p 497 
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is recommended, where is the friction factor as produced by stationary, growing 

drops on the wall. As can be surmised, is a function of and the vapor Reynold's 

number ft . This relationship is established in Reference 48 and is represented in 

Figure 26. is determined by means of Equation 35, This approach for accounting 

for the combined wall shear and drag of drops on the wall is based on the following factors: 

1 . Mercury drop populations, estimated on the basis of heat transfer considerations, 

^ are such that approximately 40% or more of the internal condensing area is 

covered . 

2. Similar area coverage has been observed in the dropwise condensation of steam 
(Ref 46). With such large area coverage, Schlichting (Ref 47) has shown that 
spheres have the effect of a rough surface and Nikuradse's data for artifically 
roughened pipes, reproduced in Figure 25 may be used. In applying Figure 25 
it is recommended that the equivalent sand roughness be assumed equal to the 
average drop size which is equal to 2/3 the maximum or critical drop size at en- 
trainment. 


The third term of Equation 32 is the pressure gradient induced by the vapor velocity changes 
caused by condensation and compressibility. For incompressible flow the momentum term 
can be calculated as follows: 
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(Eq. 39) 


Equations 33, 38 and 39 can thus be used to predict the pressure drop during the non-wetting 
condensation of mercury in zero gravity. These equations are not very usable, however, and 
a computer simplification is required to make them so. Section 7 contains a comparison 
between a simplification of the theory and some data obtained. 

The theory above gives a droplet size and number distribution function which could be used 
to calculate inventory. Consideration of the liquid on the wall is also required, of course. 

A major portion of the inventory will be in the low quality end, however, which is not 
amendable to analysis due to the influence of gravity on agglomeration. 

Pressure fluctuations at the interface of a non-wetting condenser are principally due to 
impingement on the interface of liquid that has been picked up by the vapor or dragged along 
the wall. Again, the theory could be used except for the influence of gravity. The theory 
provides the droplet mass, velocity, and number from which momentum could be determined. 
The number of drops is so great, however, that the pressure pulsations would be so frequent 
as to be undetectable . 
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5.0 RIG DESIGN, INSTRUMENTATION, AND TEST PROCEDURE 


5.1 Rig Design 

A closed, natural circulation test loop was chosen to provide the wetting and non-wetting 
experimental data required. Figure 27 is a schematic representation of the test loop and 
Figures 28 and 29 are photographs of the rig. The use of a closed loop allowed continuous 
operation, although loop dynamics problems occurred as a result. Non-wetting data have been 
successfully obtained with the rig described below. Wetting data will be obtained with essen- 
tially the same rig, with wetting condensation induced by previously conditioning the tube. 

Basically, the rig consisted of a mercury pot boiler, a heated and insulated line superheater, 
a horizontal condensing test section, a condensed-liquid collecting and flow metering system, 
and a return line to the boiler. The boiler was immersed in an electrically heated salt bath. 
Pressure taps were placed at 14.5 inch intervals on the replaceable test section for most of 
the testing and at 18 inch intervals for some of the earlier tests. The pressure measuring 
system consisted of mercury liquid-filled "U" tube manometers constructed from stainless 
steel and transparent plastic. All manometers were connected to an adjustable argon supply 
pressure to control the liquid mercury level. A shop air supply system provided a means for 
variable cooling of the condensing system. 

Argon was used as a general cover gas for the system. An x-ray generator, control, fluoroscopic 
screen and appropriate shielding were mounted on a dolly so that the full 7.5 to 8 foot con- 
densing length could be inspected and the condensing process thereby observed. This arrange- 
ment is illustrated in Figure 30. 

The possibility of air leakage into the system was minimized by using welded bel lows-sealed 
valves, argon cover gas(LindeHi Purity 99.995%) and a mercury head. These features are 
particularly important to insure that the wetting condition is maintained, as de-wetting 
between mercury and metal surfaces has been observed after exposure to air (Ref 49). 
Maintenance of a relatively air-free system is also important to maintain a non-wetted 
condition. "Pseudo -wet ting” was observed in experiments at TRW as a result of con- 
tamination due to air in the system. See Section 6.2. 

5.1.1 Boiler 

The boiler was an 8 inch diameter by 9.75 inch high stainless steel bucket partially immersed 
in a heat treat salt composed of 52% potassium nitrate, 40% sodium nitrite and 8% sodium 
nitrate. This salt melts at about 400°F and can be used in the range 500 to 1000°F. Electric 
immersion heaters capable of 9 kilowatt power output were used to heat the salt bath and 
thereby control the mercury flow rate. At mercury heights of 6, 5, and 4 inches, about 140, 

118, and 93 pounds of mercury respectively, were contained in the boiler (at 700°F, = 

795 pounds per cubic foot). The salt level was maintained below that of the mercury but such 
that 3 pounds per minute of mercury could be evaporated at a heat flux of less than 20,000 
Btu/per hour via the heat transfer area. 
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X-RAY CONTROL BOX 
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FIGURE 27 


MERCURY CONDENSING TEST RIG SCHEMATIC 






MERCURY CONDENSING TEST RIG (BOILER END) 
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FIGURE 28 




MERCURY CONDENSING RIG (FLOW METERING END) 
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FIGURE 29 



MERCURY CONDENSING RIG (X-RAY GENERATOR) 
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FIGURE 30 


Unfortunately, flow limitations coupled with excessive salt temperatures and pressure oscillations 
were encountered. The flow limitations occurred at low heat transfer areas and low mercury 
levels and the pressure oscillations occurred at high mercury levels. Careful adjustment of 
the salt and mercury levels was required to avoid these problems. At heat fluxes approaching 
and greater than 30,000 Btu/hr-ft^, salt temperatures were unsteady and sometimes would 
unpredictably start to Increase. Possibly partial film boiling occurred either on the total 
or on only a portion of the boiling surface. Liquid level in the boiler was maintained at 
four Inches or less to minimize the possibility of a head causing a sufficient variation In 
saturation temperature to result in different boiling mechanisms (l.e., nucleate, partial 
film and film) occurring at various portions of the boiling surface. 

A baffle-cup arrangement was provided at the vapor exit of the mercury pot boiler, as shown 
in Figure 27, to minimize carry-over. A drain iine from the baffle cup was provided for the 
return of liquid that dropped from the vapor. The vapor travelled from the boiler to the 
condensing section entrance via one-inch outside diameter by 0.035 inch wall superheater 
tube. Vapor velocity at the boiler exit was thus less than 20 feet per second. Several 
previous experiments conducted at TRW with transparent boilers and condensing sections 
showed no visible evidence of carry-over from pot-type boilers at this velocity. 

5.1.2 Superheater-Entrance Section 

The one-inch stainless steel tube exiting from the boiler rose about 4 feet vertically and was 
brought to a horizontal position by a 6 inch radius bend. A reducer was employed in the 
transition from the one-inch tube to the condensing section entrance length which is approxi- 
mately one foot long. The inside diameter of the entrance section was matched with that of 
the condensing test section. Welding couplings were used on either side of the entrance 
length to avoid any weld penetration and subsequent flow interference. 

The complete superheater-entrance section was heated and insulated to prevent condensation 
before the condensing section was entered. The heaters were electrical resistance Chromolox, 
KTL 1-390 WA, 4500 Watts, 120 V, 1 phase, circular immersion heater with 24 inch riser 
type with a stainless steel armoured protective covering. An immersion thermocouple was 
installed just below the bend in the one-inch tube to measure the degree of superheat. 

5.1.3 Condensing Section 

The condensing section geometries were described in Section 3. All test sections were 
horizontal within 0.5 inch over the 8 foot length (within 0.3 degrees of angle). The test 
section was kept straight by 12 supports tied directly to the air manifold. Each support was 
fixed to the test section by three screws to give nearly point contact and thereby minimize 
heat loss to the support. Previous to installation, each section was treated as follows: 

a. cleaned with trichlorethylene, 

b. steam cleaned, 

c. hot air dryed. 
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The 1/8 Inch by 0.020 Inch wall stainless steel manometer taps were then welded to the bottom 
of the condensing section. The spacing of the local pressure manometers, from the beginning 
of cooling, was follows: 

Test Series A Test Series D, E, and F 


No. 

Location, inches 

No. 

Location, inches 

1 (and A) 

0 

1 

0 

2 

18 

A 

0.5 

3 

36 

2 

14.5 

4 

54 

3 

29.0 

5 

72 

4 

43.5 

6 

90 

5 

58.0 

7 

102 

6 

72.5 



7 

87.0 



8 

101.5 


NOTE: "A" denotes absolute pressure measurement 

In welding the manometers to the condensing section, the amount of filler rod was minimized 
to reduce heat loss down the manometer taps. This was desirable in order that the normal 
mercury condensing pattern would not be disturbed by cold spots in the condensing section 
at the pressure tops. 

To obtain system pressure level, an absolute pressure reading was obtained from the first 
pressure tap for the Series A tests. This reading was also used to measure local pressures. 

For the D, E, and F series a separate tap was provided for this pressure measurement as the 
series A testing showed that it was difficult to keep the first tap filled when the absolute 
pressure was also read from this tap. For test series E and F a cooling water line was brazed 
about 3/4 inch down the pressure tap to insure that the taps were filled with mercury. 

For the F series, however, the cooling water line was not brazed to the absolute manometer 
tap to minimize the localized cooling that is thus produced. 

The test section was cooled by two approximately diametrically opposed air streams as shown 
In Section A of Figure 27. Cooling load was adjusted by varying the shop air supply between 
50 to 100 psig and by adjusting the air manifold spacing from the condensing section from 
1/2 to 4 Inches. The heat flux available was 25,000 Btu/per hour per square foot or less, 
which is typical for present day direct radiator-condensers for mercury Ranklne cycle space 
power conversion systems. 

Each manifold consisted of 164 1/16 Inch holes at a spacing of 3/4 inch. Parallel flow 
guiding strips 3/4 inch long and about 1/4 inch apart were provided along the full length of 
the manifold. Air entered at the top of each manifold at four locations spaced to provide 
an even distribution of air flow. The air entered the baffle via a partially inserted closed 
end tube with holes around the side. This avoided an increased flow through the jets In the 
immediate vicinity of the manifold air entrance. 
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5.1.4 Fluoroscopic Observation 


An x-ray* -fluoroscopic screen combination mounted on a dolly was used to observe (1) 
some of the internal phenomena occurring during condensing, (2) the liquid level in the 
manometers, and (3) the vapor-to-liquid interface. Transportation of the assembly was 
possible over the full length of the condensing section. The x-rays passed through the 
condensing section to a 4 by 4 inch fluoroscopic sheet. 

For safety a lead sheet covered the personnel side of the x-ray generator. In addition, the 
fluoroscopic sheet was mounted behind two sheets of 1/4 inch leaded glass and centrally 
mounted in a 2 foot by 2 foot by 1/2 inch lead sheet. Personnel also maintained a safe 
distance when the x-ray generator was operating. Victoreen, 0 to 200 mill! -roentgen 
"Direct Reading Dosimeters" were carried by ail personnel to monitor radiation pickup. 

The resolution of the x-ray secreen observations was such that low-velocity mercury drops of 
about 0.020 inch were discernable within 0.028 inch stainless steel and 0.010 Haynes 25 
tubes. Much smaller drops could be observed if an x-ray fluoroscopic combination with 
better resolution were employed. 

5.1.5 Condensing Length Control and Return 

Completion of condensing was always taken to be the point where the tube was observed on 
the fluoroscope to be completely full. Vapor-to-liquid interface or condensing length 
control was achieved by means of a high temperature welded bellows sealed needle valve 
In the subcooling line between the interface and the flow metering sight glass. Additional 
control was possible by varying the back pressure on this valve via an argon pressure source 
and vacuum pump combination. A plenum chamber was provided to which argon could be 
added or removed to maintain a pressure differential across the valve to provide controllability. 
The reservoir could also be used to store excess system Inventory if necessary. The possible 
effects of argon (as a non -condensible) were reduced to insignificance by the rig operating 
and data taking procedures; no appreciable quantities of gas were observed via the 
fluoroscope, and the thermocouples inches before the interface were periodically checked to 
be near saturation temperature. 

The mercury return line included a provision for system draining and Installation of a pump 
to circulate the mercury with additives to promote wetting (see Section A, Appendix A). 


5.2 Instrumentation and Data Accuracy 


*Picker Corporation, Cleveland, Ohio, Portable X-Ray Model 6170, 110 KV Hot Shot with 
40-110 KV stepless control; 1-10 ma stepless control; 1 10 V, 1 phase, 60 cycle power; con- 
tinuous service at 3.5 ma and 110 KV; 0.5 mm focal spot; no backward radiation. 
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5.2.1 Pressure Measurement 


Mercury-filled U-tube manometers were used to measure pressure at the condensing section 
locations listed in Section 5.1 .3. As a result one absolute and several relative pressures were 
provided. Figure 27 schematically depicted these pressure taps while Figures 28, 29 and 
30 are photographs of the rig and show the manometers installed. 

Each leg of the U-tube manometer consisted of a 1/8 inch by 0.020 inch wall line with a 
welded bellows sealed valve about one foot from the test section. Closing of these valves 
allowed variation of system pressure without mercury spill-over. After a bend in the 1/8 
inch tube, transparent plastic tubing (1/4 inch outside diameter clear plexiglass) was coupled 
in to allow the reading of liquid level on meter sticks. The plastic tubing (1/4 inch inside 
diameter "Tygon" vacuum tubing) was then continued to a manometer pressure control 
chamber connected to an argon supply and vacuum pump. Variation of the manometer back 
pressure maintained the liquid level such that it could be read on the meter stick. 

The liquid level in the manometer near the condensing section was maintained by either 
liquid falling down into or vapor condensing in the manometer. Observation of the liquid 
level on the fluoroscopic screen indicated that at a relatively steady pressure level the liquid 
level could be 1/2 inch below the test section. For test series A and D, then, the manometer 
control pressure was always increased to refill the lines. Pressure levels were then read after 
about 30 seconds to allow the system to return to equilibrium. In this time the liquid level 
remained within about 1/32 inch of the test section. 

For test series E and F, a water cooling line was introduced about 3/4 inch down the 
manometer. Subsequent fluoroscopic observations showed manometer fillage to be within 
1/32 inch (0.9 millimeters of mercury) of the condensing section. Allowing for some error 
in reading the liquid manometer and in the relative alignment of the meter sticks, an 
accuracy of * 3 millimeters of mercury in reading each local pressure is estimated if system 
pressure is stable. Pressure oscillations of the order of ± 5 millimeters were observed. These 
have been noted appropriately. The best accuracy for any pressure drop measurement would 
then be * 6 millimeters or less than * 0. 12 pound per square inch plus an error due to 
possible tilt in the tube of * 1 .59 millimeters per foot between the pressure taps. For a 
14-lnch increment then, accuracy was ^0.15 pound per square Inch plus what may be 
noted on the data tables in Section 6.0. 

To eliminate gases from the U-tube manometer, a vacuum was kept on the system for at least 
two hours before startup and the liquid in the metal leg of the "U" was replaced by mercury 
from the transparent leg which was first observed to be gas free. Errors due to gas entrain- 
ment in the manometers were thus prevented. 

5.2.2 Flow Rate 

Flow rate was measured by collecting and timing mercury exiting from the condensing test 
section in a graduated flow metering sight glass, as shown in Figure 31. Closure of the valve 
immediately downstream caused liquid buildup which was timed. The liquid buildup varied 
from to 0.84 to 1 .69 pounds of mercury. Time intervals varied from 40 to 90 seconds. 
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During the flow measurement, the interface control valve usually required modulation to avoid 
any change in subcooler inventory, thus leading to an erroneous flow rate. Interface position 
was observed on the fluoroscope and held to within about +1/4 inch. As a result/ an accuracy 
of about + 0.6% in flow could be introduced. The accuracy of the flow was also affected 
by the location in the sight glass at which time was started . At large mass increments, time 
was started in the very narrow lower portion of the flow metering sight glass, thereby re- 
ducing the mass increment error. The estimated accuracy for the various flow increments 
obtained and presented on the data sheets of Section 6 is indicated below as a function of 
mass increment and interface diameter: 


Mass Increment, pounds Interface Diameter, inches 



0.397 

0.319 

0.2 

1.67 

1.1% 

0.8% 

0.4% 

1.54 

2.1 

1.8 

1.4 

0.84 

3.3 

3.0 

2.6 


1 

Assumptions: 



(1) Interface varied + 1/4 inch between 
start and finish of measurement. 

(2) Liquid Level related to time to within 
1/16 inch. 

The accuracies are based on the assumption that interface was held between +1/4 inch and 
liquid level could be observed and related to stopwatch time to within + 1/16 inch. Both 
these values are pessimistic in the sense that better accuracy than that presented was very 
often obtained. Whenever possible the data were obtained at the largest mass increment 
(which was "99 and 44/100%" of the time), thus minimizing the error in flow rate. 

5.2.3 Temperature 

Temperature measurements were required in the following key areas of the condensing rig: 

a) One-inch upstream of the interface to provide a check for non-condensibles. 
(Several condensing experiments at TRW have measured tube wall temperatures 
hundreds of degrees below saturation due to non-condensibles.) 

b) salt bath temperatures monitored for overheating and steady state, 

c) superheat. 
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All temperatures where measured with ^24 Duplex iron constantan thermocouple wires. The 
thermocouples were spot welded to the area of interest except for the superheat immersion 
probe. Superheat temperature was read out on a Thermoelectric Corporation ‘'Minimite" 
with 2°F intervals over a 0 to 1000°F or 900 to 1800°F range. All other temperatures were 
read out on a temperature logger. 

5.2.4 Vacuum Equipment 

No vacuum gage as such was employed in the system. Vacuum levels are of importance 
when dealing with surface chemistry phenomena such as wetting and non-wetting. The 
vacuum pump used in the system is then worthy of description. A WM Welch Scientific Co. 

Duo Seal 1405 B, 0.1 micron vacuum pump was used. 

5.3 Rig Operating Procedures 

The procedures for startup, approaching the data point, taking point, and both emergency 
and normal shutdown are presented below. These procedures have been employed during 
the non-wetting testing and are expected to be applicable to the wetting tests. For aid in 
understanding these procedures reference to Section 5.1 is required. 

To start testing, all manometer valves were closed and the boiler heaters, the argon supply, 
and the vacuum pump were turned on. After the system had been evacuated, liquid in the 
1/8 inch stainless steel tube of the manometers was replaced by opening the manometer 
valves. All manometer valves were then reclosed. As flow began, all airmanifold supports 
were inspected for twisting caused by thermal expansion to prevent errors in pressure readings 
due to gases in the liquid legs. The system was allowed to run for 15 minutes with the flow 
control valve open to remove gas. The cooling air supply was then turned on and the vacuum 
valve to the plenum was closed. The interface control valve was then closed to allow the 
test section to fill with mercury condensate. Before approaching the data point, all manometers 
were refilled with mercury condensate from the test section. 

After the desired data point conditions had been determined, T $a | |. - T^gsat was estimated 
from past experience for the desired flow. This AT was then added to the desired saturation 
temperature to obtain an estimated boiler bath temperature. The boiler temperature was 
adjusted to the desired level as quickly as possible by manipulating air flow and boiler 
heaters. When the boiler heater was within 10 degrees of desired temperature, the boiler 
heaters were set to produce the desired mercury flow rate and the air manifold setting and/or 
air flow rate were adjusted until the desired condenser pressure was reached. Condensing 
length was controlled by means of the interface control valve and back pressure. The 
mercury flow rate was periodically checked as the system stabilized and manometers were 
opened when stability was attained. 

Before taking data, the following stability criteria were met: 

1 . The interface was under control near the final position for at least 15 minutes. 
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2. Logger boiler bath temperatures were steady (less than 3 degrees change in 15 
minutes). 

3. Air manifold pressure remained steady. (Pressure regulators were eventually used.) 

4. Flow rate was within 5% of the desired flow. 

5. No gas was observed at the interface on the fluoroscope screen and measured 
temperature before the interface read ~ T sa j.. 

6. System pressure held within 20 millimeters of the final pressure for 15 minutes 
and was within 10% of that desired. 

Liquid level was checked to insure that the manometers were full. Manometers were read 
only when full and holding steady, never during a time when system pressure was Increasing 
and manometers could not be filled to the top. Immediately after the pressure readings, 
the sight glass valve of the collector was closed and flow rate was read. The timer for flow 
rate readings was started as the mercury level passed the top of the lower reference band, 
allowing for buildup to the upper reference band. Careful maintenance of the interface 
posttlon in the condenser was required during this time. The timer was stopped when the 
mercury level reached the top of the upper reference band. Bands used and buildup time 
were recorded. 

After at least two data recordings had been made at each point and checked for similarity, 
preparations were made to change to the next data point. Manometer levels were carefully 
maintained as system pressure changed. 

To terminate the test, all heaters were shut off, the flow control and return line valves were 
shut, and the test section was allowed to fill with liquid condensate. To prevent excessive 
pressure buildup in the boiler, the boiler bath temperature was reduced to below 800°F 
before shut down. All manometers were then refilled to 800 mm, the x-ray was shut off, 
the vacuum valve on the collection pot was opened, and the closed-off section of the system 
was then evacuated. 

The procedure established for emergency shut down consisted of four steps: 

1 . Shut off all heaters and the x-ray as quickly as possible. 

2 . Shut manometer val ves . 

3. Shut flow control and return line valves but leave the air manifolds on. 

4. Take any possible steps to correct the situation. 

All safety procedures that have been established for mercury, x-ray, and fire hazards were 
observed throughout the tests. 
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6.0 EXPERIMENTAL RESULTS 


6 . 1 Non-Wetting Condensing Data 

Experimental non-wetting data are presented in Tables II, III, IV and V. The absolute 
pressure at each local position is presented . Pressure tap location is identified by distance 
from the beginning of the air cooling which is where the test section inlet insulation and 
guard heaters end. For example, P 54 is the pressure 54 inches from the beginning of con- 
densation. Actually the absolute pressure for the zero distance only was measured and all 
others were measured relative to this location and to each other, as described in Section 
5 . 2 . 1 . Because relative pressures were measured, a correction was required across the 
vapor-fo-!iquid interface only. Pressures measured by the first pressure tap in the liquid 
portion of the tube have been corrected by substracting a head of one-half the test section 
tube diameter. For the 0.319 inch data then, a correction of -4 millimeters was required. 

The accuracy of each pressure measurement is discussed In Section 5.2.1. System oscillattons 
caused pressure fluctuations which were recorded and are presented in the data tables. 

Flow rate was calculated from the mass increment and collection times presented. Flow 
accuracy as a function of mass increment collected and tube diameter is estimated in Table 
II. For the Series A, D, and E data, two corresponding sets of pressure and flow data were 
obtained at approximately the same system pressure level and flow rate (e.g., A-3 and A-4). 
For the F series tests the two sets of pressure profiles were recorded and the flows were 
averaged and associated with each pressure profile. 

To give an indication of the superheat existing, the me asured inlet temperature is presented 
along with the saturation temperature (T $a j.) corresponding to the test section inlet pressure 
(P Q ). Superheat is then (T sa f - T 0 ). Condensing length (L c ) , the approximate axial length 
of the interface (L|) , and the distance from inlet (L') at which globules were fluoroscopically 
observed on the tube bottom are also given. L| and L' were not obtained for the A, D, and 
some of the E series data. 

Table II contains some data (W series) obtained during the "pseudo-wetting" experienced 
with the stainless steel 0.319 inch tube, described in Section 6.2. 

Figures 32, 33 and 34 summarize the range of flows and saturation levels covered by the data. 
Details pertinent to the apparatus, instrumentation, and operating and data taking procedures 
were presented in Section 5.3. 

6.2 1 'Pseudo-Wetting 11 Investigations 

During initial testing of the 0.319 inch inside diameter condenser tube (A series), partial 
wetting was observed in the low quality region of the tube after about 70 hours of operation 
tion. The observations were made on the fluoroscope screen. Figure 35 shows sketches of 
the appearance of wetting and non-wetting condensing on the fluoroscope screen. 

The argon cover gas was replaced by an air cover in an attempt to obtain more data on the 
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TABLE I 

MEASURED DATA DURING NON-WETTING CONDENSATION OF 

SERIES A AND W (0.319 inc 


Hg Collected, 
lb 


Ib/rnin. I mm Hg 














TABLE III 

MEASURED DATA DURING NON-WETTING CONDENSING OF MERCURY IN A HAYNES 25 TUBE 
SERIES F (0.319 Inch ID by 8 Feet Long) 
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♦Approximate oscillation of pressure readings. 

*Not measured but obtained from Pq and saturated -temperature relationship. 

Note - Subscript of P's denotes distance of manometer tap from the beginning of cooling. 
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MEASURED DATA DURING NON-WETTING CONDENSING OF MERCURY IN A HAYNES 25 TUBE 
SERIES E (0.397 Inch ID by 8 Feet Long) 
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Approximate accuracy of pressure readings. 













FLOW, PRESSURE AND CONDENSING LENGTH DATA POINT COMBINATIONS 
OBTAINED FOR THE HORIZONTAL CONDENSATION OF MERCURY IN A HAYNES-25 
SPINEL COATED 0.4 TO 0.2 INCH ID BY 7 FOOT LONG TAPERED TUBE 

(SERIES D) 



FIGURE 33 






FLOW, PRESSURE AND CONDENSING LENGTH DATA POINT COMBINATIONS 
OBTAINED FOR THE HORIZONTAL CONDENSATION OF MERCURY IN A HAYNES-25 
AND A STAINLESS STEEL 0.319 INCH BY 8 FOOT LONG TUBE 

(SERIES A & F) 



FIGURE 34 
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non-wetting condensing of mercury. Further testing, however, showed that wetting still 
existed and progressed with time. Since the existence and development of the wetting is 
of basic interest, observations were made and data taken for several days before termination 
of this test. The above performance will be compared with the performance of tubes with 
wetting condensing promoted by additives. Figure 36 compares pressure profiles during the 
wetting and non-wetting conditions and compares each with Martinelli theory (See Section 
4.2). 

Continuing with the non-wetting testing, a 0.397 inch inside diameter by 8 foot Type 304 
stainless steel condensing tube was installed in the rig. An air cover gas was maintained in. 
order to maximize the non-wetting testing time. Contrary to our expectations, wetting was 
noted within several hours. Similar results were obtained with the 0.4 inch to 0.2 inch 
inside diameter by 7 foot Haynes-25 tapered tube. These results are very unusual when one 
considers that de-wetting between mercury and metal surfaces has been observed after ex- 
posure to air (Ref 49). Introduction of air into the loop was expected to create this type 
of contamination and thereby cause the mercury to de-wet the surface. 

The premature wetting obtained can be attributed to 'jsseudo -wetting" between mercury 
saturated contaminants or crud and contaminant on the surface. The problem was to determine 
the cause of the "pseudo ^wetting" and how it could be prevented. In order to continue 
the testing and insure that sufficient non-wetting testing time would be available, the following 
specific areas were investigated: 

1 . the causes of the premature (or "pseudo-wetting") wetting, 

2. methods for eliminating premature wetting once obtained, 

3. steps required to insure that non-wetting data could be obtained continuously 
for at least a period of two weeks. 

To expedite this investigation three minature versions (Minirigs) of the full scale test rig 
were made. A sketch of the rig is presented below. 
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Disd 


COMPARISON OF MEASURED VARIATION OF PRESSURE WITH LENGTH 
WITH MARTINELLI FOR THE COMPLETE WETTING AND NON-WETTING 
CONDENSATION OF MERCURY IN A 0,319 INCH ID, 316 STAINLESS 

STEEL HORIZONTAL TUBE 

• Run A-26; w = 2.08 Ib/min (No wetting observed) 

■ Run W-4; w - 2.1 1 Ib/min (Wetting* observed) 

wmmma Martinet I i: See Appendix C 
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FIGURE 36 
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Operation of the rigs was fluoroscopically observed to detect any significant changes in 
wetting in the condensing section. Table VI summarizes the operation of the rigs and compares 
the operation with the full-scale rig. 

The three Minirigs were run a total of about 1000 hours. All three rigs were purposely dirtied/ 
then cleaned, and the condensing section was replaced and rerun to determine if non-wetting 
operation were possible after the rig had been contaminated. Minirig No. 1 was contaminated 
by Introduction of air after operation for about 500 hours. No wetting was observed on the 
fluoroscope screen during the 500 hours operation whereas the introduction of the air caused some 
wetting. A cleaned but otherwise untreated Haynes-25 type material was used in the condensing 
section after the rig was cleaned. Non-wetting operation was achieved for about 150 hours. 

Minirig No. 2 was contaminated by using dirty mercury and large amounts of air. The rig 
was cleaned and a Haynes-25 type material spinel-coated tube was installed in the condensing 
section. A green spinel coating is formed by high temperature ( **+ 1750°F) exposure in a 
wet hydrogen atmosphere. Non-wetting operation was obtained for about 100 hours. The 
rig was shut down and argon back filled over night and during weekends. Minirig No. 3 
was also dirtied with dirty mercury and air. The rig was cleaned and a stainless steel spinel - 
coated tube was installed in the condensing section. Due to poor cleaning the spinel coating 
was not uniform, however, and possibly caused the wetting that was noticed upon initial 
startup. A Haynes-25 type material spinel -coated tube was then used in the condensing 
section. Non-wetting operation was obtained for about 160 hours. 

As a result of the Minirig testing the following procedures are now required to prevent 
pseudo -wetting while obtaining non-wetting mercury condensing data. 

1 . Air-mercury exposure in the rigs must be minimized. 

2. When necessary, rig components should be cleaned with a dilute (20%) nitric 
acid solution followed by water. 

3. A Haynes-25 type material tube treated in a wet hydrogen atmosphere should 
be used to provide non-wetting condensation for about 100 hours if (1) and 
(2) are heeded. 

4. Rig shut down should be minimized but if necessary a high purity argon cover 
gas should be used . 

6.3 Reduced Data 


A digital computer program has been developed for the local and over-all reduction of 
constant tube diameter data. Data inputs are total condensing length, flow rate, absolute 
pressure level, local pressure readings and geometry. Average local and over-all saturation 
pressures are then calculated and thermodynamic properties are calculated according to 
Reference 50. 
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Local and over -all measured two-phase pressure drop ( AP tp ) is calculated, corrected 
for vapor momentum recovery to give total two-phase pressure drop (A P Tpr ), and divided 
by the vapor-only pressure drop ( A P ) for the same geometry and temperature. Vapor 
velocity is assumed to vary linearly with length in the data reduction. The ratio 
APtpt/* p <s can then be used for experimental correlation and comparison with 

Martinelli (see Section 4.1 .2). In addition, local liquid and vapor Reynold's numbers are 
calculated and the appropriate term calculated according to the equations of Section 
4.1 .2. An additional term Re 0 ( v jg/jo) , ^ S is calculated as this combination has 
successfully correlated over-all two-phase pressure data at TRW (Ref 3, 5). 


The computer program cannot presently handle the tapered tube data. Additional correlating 
parameters calculatable from the inputs can be easily included. The computer program is 
obtainable from either NASA or TRW. 
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PSEUDO-WETTING OPERATING HISTORIES 
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PSEUDO-WETTING OPERATING HISTORIES (Continued) 
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7.0 COMPARISON OF DATA WITH THEORY 


The theory available for comparison with mercury condensing data Is presented in Section 
4, "Theory and Analysis." At this time comparison between the data and theory is in- 
complete. The comparisons available are nevertheless presented, although a more exten- 
sive comparison is to be expected in the final report. 

Section 4.2 presents a method for predicting pressure drop during the non-wetting condensing 
of mercury with low Froude numbers (i.e., when the ratio of gravity to inertia forces is 
small). The most difficult phenomenon to account for in the non-wetting pressure drop pre- 
diction is the effect of drops in transit. A rather complex computer solution is required 
for the complete solution. To provide a check on the theory, the equations were made 
amenable to hand calculation by means of several simplifying assumptions. Figures 37 and 
38 compare the pressure profiles actually obtained with the profiles predicted by the theory 
of Section 4.2 simplified by means of the following assumptions: 

1 . Droplet drag can be determined from the Stokes equation, 

r . 2A _ 24- xty , 

4 " ' ix & T |i-e| 

2. Slip can be assumed to vary linearly with droplet size from zero to nearly one, or 

e - ^ ~ ^ — . 

■Se" 

3. The effect of the axial vapor velocity variation with length was neglected and 
Equation 37 became a linear differential equation. 

4. Only drops having slips from £ = 0 to € = 0.96 were considered. 

The limit of € = 0.96 is effectively a correlating adjustment. 


The low pressure gradient at the inlet of the 0.397 inch inside diameter tube data of 
Figure 38 is tentatively attributed to the increased cooling effect caused by the two 
water-cooled pressure taps used at the inlet during these tests. Otherwise the theory 
appears to correlate quite well with the data. 

Section 4.1 describes the application of Martinelli's data to the prediction of condensing 
pressure drop. A comparison with wetting and non-wetting data was felt to be useful; 
this is presented in Figure 36. Both non-wetting and pseudo-wetting pressure data are 
presented for the 0.319 inch inside diameter tube. Appendix C presents the sample 
calculations made to determine the pressure curve shown in Figure 36. 
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COMPARISON OF SIMPLIFIED THEORY WITH DATA FOR THE 
NON-WETTING CONDENSING OF MERCURY 

Run - A-26 O - Measured data 

W = 2.08 Ib/mm 

ID = 0.319 inch Material - Stainless steel 

Horizontal 



FIGURE 37 
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COMPARISON OF SIMPLIFIED THEORY WITH DATA FOR THE 
NON-WETTING CONDENSING OF MERCURY 


Run - E— 1 5 Material - Haynes 25 

W - 2.38 Ib/min Horizontal 

ID - 0.397 inch 



FIGURE 38 
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APPENDIX A 


WETTING AND NON-WETTING CONDENSING 

The general consideration of interfaces between solids and either liquid and/or vapor is not 
simple. Considerations are required of surface roughness, non-homogeneity, and con- 
tamination. Contamination may include very lightly held, physically absorbed molecules 
or thin oxide layers. These phenomena are summarized by Bikerman in Reference 52. 

The case of non-wetting mercury can be attributed to surface contamination (oxides) and/or 
adsorbed gases. An example of the effect of adsorbed gases and/or surface contaminations 
is the fact that mercury has been shown to wet glass that was heated under a vacuum 
(Ref 53). Similarly metal that is normally non-wetting, such as stainless steel, when broken 
in mercury to give immediate exposure to the mercury is immediately wetted. 

With sufficient time and temperature mercury will dissolve adsorbed gases and remove the 
oxide surface coating. In the interim, however, the surface will be non-wetting and non- 
wetting condensation will occur. Non-wetting data were thus obtained from tubes that were 
cleaned according to procedures commonly used for cleaning mercury Rankine cycle space 
power systems. These procedures include sufficient exposure to oxygen to allow some form of 
oxygen contamination. The oxygen pressure required to produce contamination is far lower 
than the best obtainable vacuum for most metals of engineering importance. 

Wetting condensation can be obtained by using magnesium and titanium additives to remove 
contaminants. (Ref 54, 55) Both magnesium and titanium are "getters" of oxygen and break 
down oxide films as well as remove adsorbed gases. As a result, a clean metal surface is 
left for the mercury to wet. Actually, however, mercury containing titanium forms a 
titanium intermetallic of some sort (iron, carbon, or nitrogen, and even mercury) which 
serves as a barrier to mercury contact with the metal , and mercury wets this intermetallic. 
(Ref 56). If the titanium intermetallic should be eroded, the clean metal surface remains 
and should cause wetting condensation. 

Wetting condensation can thus be maintained as long as oxygen and other contaminants are 
kept from the surfaces. To prevent reversion to non-wetting, the surfaces to be preserved 
should be exposed to mercury vapor and/or liquid only. As a poor alternative an inert gas 
cover such as very high purity argon might be used. 

Another form of wetting has been obtained in TRW experiments and can readily be observed 
by placing contaminated mercury in a glass cup. A non-uniform and spotty mercury-like 
substance will adhere to the glass. This form of wetting is pseudo-wetting in that it is really 
wetting between mercury-soaked contaminants and the contaminated glass wall. Section 6.2 
discusses the actual occurrence of "pseudo-wetting" during test. An obvious problem thus 
exists. When "wetting" appears. Is it actual wetting or not and how do wetting and pseudo 
wetting compare. This is a surface chemistry problem beyond the scope of this contract 
but certainly worthy of consideration. 
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APPENDIX B 

NON-WETTING DROP MODEL DERIVATIONS 
DROP DISTRIBUTION FUNCTION 


Flow Rate 


Let be the number of drops per unit volume, having diameters between & and 

d& , at distance l_ . 

s' 

The mass of one drop of diameter 6 is . 

O 

The mass of fluid per unit voulme, contained In drops having diameters between 6 and 
J6 at distance L , is 


The mass flow rate, at distance L. , of fluid contained in drops having diameters between 
d and & is then 




Let U + (i N 0=f(i,L)U v (L^ 


Equation B-i gives, 

U v (L)N(^LH(<S A 0^dS 

The total mass flow rate of fluid contained in drops at distance L is 


u»0-> f N«,L) f 

It is known that d ^(Cj/d L is always positive. 
c = ^ 


Let 


(Eq B-1) 
(Eq B-2) 

(Eq B-3) 
(Eq B-4) 

(Eq B-5) 
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(Eq B-6) 


Thus, 


UJ 4 ( L'l r C,u v (L)/N(4U€(4,U)& 3 d<5 




Sketch Showing Variation of Drop Size With Length 


jQl 




o 

o 




L=o 

4(0=4, 


L V-4<lL 

4 to 4ti-+<*0 = 

& c co 4* d^-CO = 

d/O -f 


a w 

The rate of entrainment of liquid by vapor is — M r 

a L. 

The increase in drop flow rate between L and L +<JL j$ 

d<*V L > = 

This Increase of flow rate consists of drops having diameter between and 

and includes all drops in this size range. 

At distance L+dL the mass flow rate of drops having diameters between <£ c and 
4+cU c is, from Equation B-3 f 

c,u v (OKiC4 e ,Ue(«s e ,L)^<IA c = 


or 


A uixC O 


duJ-f\ 


Since ,C, ^ n U v 0-') 

d L 


C, U v ^4^j. L > C, u v ( L) J <4 / (Eq B-7) 


and are all positive finite 

a l— 


quantities the product N €. (<4^ *-) must be positive and finite. But it is 

known that € therefore the values of the distribution function M (4 n l^ must 

approach infinity as the drop size approaches the critical size for that particular distance L. . 
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Continuity Relation 

The mass rate of flow of drops in the size range i to ^ + d<J is the same at all values 
of l- provided that this size is less than the critical size in the range of L considered. 
In equation form, 

dlU* ” d UJ + L £ c (Lj\J 

or, * 

C x ujL)H(S^i t )^(6 s L,)S d£ = C,u w (ON(<S,Q66S n Lj&<li 

Therefore, 


u v (L ! > >KK£,Lt)<:(<J,L^ = Uv(l 2 )K|(<S^ L^f ^ 


At L = L c , l.e. where d® £ c (L) f from Equation B-7, 


N(4 e )^^0- 


C, UyO-cU’ 



(Eq B-9) 


N(^L)fe^L) 


J_ _} 

C,S 3 U. y (0 


d 

dL 

d<k(Q 

dL 


L=L C 


(Eq B-10) 


NOTE: The quantity constant for all values of L. 

value of <5 .It is not constant for all values of 6 at a fixed vlaue of L 


Jl 


4 d, CL > 

d L 


is constant for all values of d 


L«L, 


at a fixed 
unless 



In the expression for O f Equation B— 10, the derivative 

is not a directly measured quantity. However, since d -d 

T 


^jCO> . dd c <ud _ _ 4- 

d L dL d Ul» "TTli^e^ 


A UXf(Q d^cCU.'i ^ 
d L d (ty 


(EqB-11) 


Therefore, in Equation B— TO, 

A tu^CQ d uj f C U 

d L _ dL _ _ Trt^/V 

~ 4 <WO ddcCay^ ~ 4 r dd<(Uv fl 

"3T” /*y d L d Uy |_ duiyf J 

Also, 

J_ TTbVv _ 6/°y 

<=, 4 - TTft. 

And, UJyC'-C) * 

t _ UvG-d _J — UJy('-fc) _l _ mp _ XCL^ ) I 

U V (LT ~ UJ V (0 U/Ld" x(l') Wy(Q 

^Jo 

Thus, Equation B-10 becomes. 


fcj/f l\- <oP* > 1 _ 

' ' -n-fif 4 s e «,0 xCO r u , d£s<i^ 

L c d <iv c J, 


(Eq 6 - 12 ) 


Alternatively, without involving the quantity ■* which is known only indirectly as 

a function of Z 


v | / jr . \ _ 6 ^V * | 1 

NW » 1 t ft a/Ll r d5c(^1 

L <U„ J 


Or, in terms of the quality at V_ J 


/ttdVv^ 

, 1 

1 

(• 4 .UU J 

1 ^eW.L'ixCOj 

dlly J 


(Eq B— 13) 


(Eq B— 14) 
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DRAG ON A SINGLE DROP (of diameter i at distance t- ) 


p,d — -o — 



+ direction (<v\creism| L) 


But, 


F sd = c d A s 4 


fv ( u r»l) - Cd A s< | F '* — Milliil ~ u d 


U v - uu r u v (I - ^)= Uv(i - fc) 

A s 4- ^ 


F t J - TX fy uj Cji (I - «)!• - 6| (Eq B— 15) 

2 . 

(l-€)h - € I is used rather than (|- t) to account for proper sign when 
Newton's second law of motion is 


rj Msd <1 L _ 



The mass of a single drop 



Combining Equations B— 15 and B— 16 yields 

d € j_ € dlAv _ 3 Pv Cd ( I ~ O h ~ ^ | 

dL + Uy^L - 4%T € 


d€ _ J5_ fa Cd 1 1 *c| j_ d LL y £ 
dL ‘4/» f T € “ UydL ^ 

* duv dx \ J_ dx 

w vdL l^xATT^dLrK dL 


(Eq B-17) 


(Eq B— 18) 


(Eq B— 19) 
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Equation B— 19 becomes, 

3 A Cd(«S,Q[l - t(j,L'Q|l-6(^U l _ dxXt-Y ,NJEq B-20) 

dL ^ ^ <5 €(^^U nCOdL 1 


Assume C«1 is a function of the Reynold's number based on drop diameter and vapor 
properties: 


Uvlt-€l& 


R « 4 = A 


feu-t 


Re* = js v x(a|i-€(4,0te 




then. 






a*, 
* ■ 

- 3 /^v Cd 

O-Oh- 

i 1 _ 1 


dL 

"4/^ ^ 

e 

X dL 

where 

X - 

xCO 



and 


Cd< l,' 

x(l) , e 

,4) 


(Eq B— 21) 


For a given condenser tube, handling a given fluid, at a given flow rate per unit area, a 
number of solutions must be obtained for various drop diameters. 

Start of Solution at 6 = 0 


There is a problem in starting a step-by-steo numerical solution of the differential equation 

(Equation B-21), because the derivative 2§ is infinite when €. is zero. 

d 


However, for small 

d€ 

dL 

This equation may be put in the form 
dfe!) - £ C do 


the equation has the form 

€ 


3 Cdc, > 

X € 


-3U 

with solution 


T 


\^(°4 


Cdo ^ 


7 ?TT 



which may be used, to start the step-by-step solution (see sketch) 
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TOTAL FORCE BETWEEN DROPS AND VAPOR (at distance L) 
The total force acting on a single drop has been shown to be. 


F S< J “ “7 Cj ^ (l "* €.^\\ - €\ 


(See Equation B— 15). 


The force per unit volume on drops having diameters between <& and S. + 4& Is then, 

d F j(<S,0 = F sd <4,L'> N (4,0 

- H pvu^(aCj(S,ON(S,C)[v W,U]|l - t(4,L)|4 di 
®I C 

The total force per unit volume on al^ dr<^ps at 1 is 

F d M = g pv4w /cj*(4,L)WW,o[i - e(<S,ij)|l -«W,l}|S*<jS 

J& 

Effect on Pressure Drop 


P+ JP = PF 
dV = AdL 


F,dV 


+ Jirodtoh 

P+4P 


If there are no effects other than droplet drag, 

F v dV = AdP 
J^AdL =A ^dL 
&r) drag ^v iw 

- X*(l yCiUXl N(4,L)[|- f (S,Q)|l-t(4,iJ|4*di 

‘ * Jr (Eq B-22) 


The equation for the pressure drop due to drag has been written in terms of the drag co- 
efficient which, in view of the non-steady relative motion, is likely to be a 

very complicated ( and perhaps unpredictable) function of local conditions. Since in an 
analytical treatment the drag equation must be solved first for every drop size, it is probably 


8? 



better to express the pressure drop due to drag In terms of the resulting drop velocity dis- 
tribution, i.e., in a form which does not explicitly involve the drag coefficient. Sub- 
stituting from Equation B-17, 


CjO-Oli-cIA - 


3 f A ^dL^* dL <V 
X 7T 


(Eq B-23) 


from the differential drag equation, ^ CO 


-^L 


UUi 


3Tr« 
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APPENDIX C 


SAMPLE CALCULATIONS FOR WETTING CONDENSING OF MERCURY 


This section applies the recommendations of Section 4.1 .2 to the prediction of the per** 
formance of a horizontal wetting mercury condenser. The approach used is based on a 
modification of the available adiabatic two-phase flow test results reviewed in Section 


4.1.2. 


No wetting mercury data as such are available for comparison. However, a comparison of 
the pressure profile resulting from the calculations below is made with some non-wetting 
mercury condensing data in Figure 36. 


To apply the modified adiabatic relations for pressure drop and liquid (or vapor) volume 
fraction, consider an incremental condensing length. 
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But, also 
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(Eq. C-3) 
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$-> R -r 'ff >5r 1 ’ *) 

The above equations ^are set up for the Martinet Si annular or mist flow regime parameters* For 
different regimes AF^ should be replaced with the two phase frictional pressure drop 
for the particular regime. 

*Thls expression is called the Blasius law. (See p 81, Ref 57) 
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Sample calculation: 

0.319 inch U^* = 3.1 cubic feet per pound 

-3 

94 inch IJj.* = 1 .25 x 10 cubic feet per pound 

1030 millimeters (Ts«r' 0 = 705°F ) 

mercury 

UJ =2.08 pounds per ^ 0 =193 feet per second ^ 

minute 

To determine the $ and R^to be substituted in Equation C-3, determination of the flow 
regimes existing in the condenser tube is required. This is accomplished by means of the 
condenser line on Figure 15 which shows the variation of flow regime performance with 
quality for the case being considered as well as another. To introduce quality, the 
following relationships were used: 


Tube inside diameter = 
Condensing iength = 

Pc = 



(Eq. C-4) 


(Eq. C-5) 


Figure 15 shows that annular flow will exist to a quality of approximately 5%, after which 
either slug or stratified flow will exist. Slug flow should exist because completion of 
condensing would fill the tube and not allow stratified flow to exist at x *0. 

Pressure drop and inventory may now be considered for the annular (quality varying from 
100 to 5%) and slug flow regime (quality varying from 5 to 0%). Pressure fluctuations for 
the slug flow regime may also be investigated. 


For the first 18 inch increment, assuming a linear heat flux with length. 


X* = 1.0 

V = 96 - 9 

* 96 


0.905 


X z = 0.809 

‘Assumed constant throughout tube. Variation in properties could be accounted for if 
necessary although this is not usually required except when pressure levels are low. 
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To define the flow regime. 


To obtain 
zero. 



4-UJt x 

Trb JX.sj 



Re v = 38,800 
= 286 viscous 

T 

Re v = 35,100 turbulent 
31 ,500 turbulent. 

and and 3C must calculated. 
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To 5 


= 1.05 x 10“ 4 

2 . 


V =0.01025. 

^ VT . 


From Figure 16, 

= 0.009 

The average liquid volume fraction is then 0.0045. To obtain , calculate 



0.00664. 


From Figure 16, 



The following values are then substituted into Equation C-3: 


CN 

• 

u 

f ei 

= 0.809 

X = 0.905 

O 

1 ! 

AL= 18 inch 

R^= 0.009 

X 

\\ 

• 

O 
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of course, 


As a result, the pressure drop is 


A 


Friction 


/ A \r = 2.02 - 0.86 = 1.16 psi = 60.1 millimeters mercury. 

— Momentum 


The chart below summarizes the remainder of the step by step calculations. The resulting 
pressure profile is shown in Figure 36. 


CALCULATION OF PRESSURE DROP 


Incremental 

Length, 

Inches 

Mean 

Quality 



X 

Flaw 

Mechanism 

$ 


AP 

vrv yr\ u 
H<* 

0-18 

0.905 

286 

35,100 

0.00664 

VT 

1.20 

0.0045 

60.1 

18 - 36 

0.719 

846 

27,800 

0.0138 

VT 

1.24 

0.0116 

42.8 

36-54 

0.531 

1410 

20,600 

0.0231 

VT 

1.29 

0.0179 

21.7 

54 - 72 

0.344 

1970 

13,300 

0.0453 

TT 

1.58 

0.0310 

-0.3 

72 - 91 

0.151 

2540 

6,050 

0.116 

TT 

1.93 

0.0757 



Average R^, for 100 to 5% quality = 0.1448 
Quality from 5 to 5% 

Assume slug flow in the last 0.4 foot of tube. Pressure drop can be calculated according to 
Equation 27. 



This pressure drop can usually be expected to be negligible (i. e . less than 1 mm of Hg) as 
it was in this case. 
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To calculate liquid inventory from Equation 29, 



Taking *- 0.025, R y - 0.82 and R 4 - 0.18. 

Pressure fluctuations can be estimated by using Figures. 20, 21 and .22. Taking X “ 0.025 
gives - 0.975. 

To obtain the maximum amplitude, assume ft - 0.9 for which X= 0.0036. From 
Figure 20, 


but 


a ^= + €r(l-xij = S.S8 ^ViSC. 


*= o.ozs 


From Figure 22 for b = 0.319 inch (8.10 millimeters), U m = (0.56 + 0.03) meter 
per second = 1 .94 feet per second. 

Therefore, for water, 

& = 0.0738 
For mercury, then, 

9k = -0.20 o —— = ~ 2.56 psi 

1.25 * 10 J * 62.4 


If 


5.58 

1.94' 


) 


1 = 0.20 pound per square inch 
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The frequency of this fluctuation, determined from Figure 21, for - 0.9, 



Y\ 


0.4 


< 5 - 58 - ’•”> 


55 cycles per second 


For the tube as a whole, the profile is given in Figure 36 . The liquid volume fraction is 


Rf = 0.95 (0.145) + 0.05 (0.18) = 0.150. 

TOT 
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NOMENCLATURE 


A Area or constant 

k Displacement of the interface 

0 Amplitude of interface displacement 

* 

0 Wave amplitude where entrainment, or plugging (E^), occurs 

1 JB /de 

e, <&/ j ©* 

£ Constant, complex wave velocity (wave celerity) 

Cj Drag coefficient 

Propagation celerity (Figure 5) 

A Differential operator 

^ Tube diameter 

6 Base for natural logarithms 

£ Constant 

Friction factor, function ta be determined 

■f 0 Moody friction factor 

p Force 

Maximum growth rate factor (see Figure 8) 

Fi? Froude number 

G* Mass velocity 

£ Gravitational constant, 32.2 ft/sec^ 

Kyf- Latent heat of vaporatization 

K Constant of integration 
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I_ Length 

L. Breakup length 

i*" 

l_ s Length for slugging 

lr\ Natural logarithm 

M Mass 

V> Pressure oscillation frequency 

N Drop distribution function (see Section 4.2, Reduced Frequency (see Figure 21) 
P Pressure 

C) Heat transfer rate 

V Radius 

R Volume fraction 

Film Reynold's number for neutral stability 
SS Stainless steel 

S Growth rate factor (see Figure 2) 

T/S Condensing test section 

LL Velocity 

U. 2 Interface velocity of liquid film 

LLo© Velocity outside of boundary layer 

U“ Specific volume 

V Volume 

V Volume flow rate 

V© Vertical velocity component of interface 

Flow rate 



^R© Reynold's number , 
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Wf Film Weber number, 

K Quality, tUv/uAy or flowing quality 

V Distance normal to film surface 


Greek Symbols 
0< Wave number, 


o< 

r 

i 

A 

£ 

© 

* 

/* 

cr 

T 

§ 

X 

V 


AIT*/* 

# • 

®/D 


Condensate loading, Ib/sec-in of tube periphery 

Drop size, film thickness 

Finite differential 

Slip u </a v 

Time 

Wave length 

Viscosity 

Density 

Surface tension 
Shear stress 


i(^/ACi 

I(ap/am 


TPP 

cmc p >->*«? 


Stream function (Equation i), slug flow friction factor multiplier in Equation 27 


Subscripts 


C 


Critical, gravitational constant 
Drag 


100 


* 

i 

F 

m-r 

l 

1 

yrv 

max 

r» 

o 

P 

r 

re 

s 

sai* 

T 

t 

TP 

T/i 

V 

v 

UJ 


Gas or vapor 
Liquid, film 

Fraction, friction 

Heater 

Imaginary 

Liquid 

Mixture 

Maximum 

Normal 

Initial , zero gravity 

Pipe 

Real 

Relative 

Static, slug flow 
Saturation 
Single drop 
Total 
Turbulent 
Two phase 

Condensing test section 

Volume 

Vapor 

Wall 
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•z Interface 

oo Great distance outside of boundary layer 

& Drop 

V Viscous 

q- Surface tension 
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